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Abstract 

A solution of the linearized Einstein and Nambu-Goto equations is con- 
structed which describes the evaporation of a certain type of rotating cos- 
mic string - the Allen-Casper-Ottewill loop - under the action of its own 
self-gravity. The solution evaporates self-similarly, and radiates away all its 
mass-energy and momentum in a finite time. Furthermore, the correspond- 
ing weak-field metric can be matched to a remnant Minkowski spacetime 
at all points on the future light cone of the final evaporation point of the 
loop. 

Short Title: Evaporation of a String Loop 
PACS numbers: 04.25.Nx, 98.80.Cq 

1. Introduction 



Cosmic strings are thin filaments of topologically-trapped Higgs field energy whose 
dynamics, in the zero-thickness limit, are controlled by the Nambu-Goto action. 
In this limit cosmic strings are effectively line singularities which move under the 
action of their own elastic tension, and radiate gravitational energy and momen- 
tum in the process. They may have played a role in the formation of large-scale 
structure in the early Universe (see pQ for a review), but their gravitational prop- 
erties are still only poorly understood, and they consequently remain problematic 



objects for general relativity. In particular, a rigorous distributional description 
of line singularities has not yet been developed within the framework of general 
relativity, nor has it been demonstrated that zero-thickness cosmic strings can al- 
ways be interpreted as the unique limits of well-behaved finite-thickness solutions 

At present, all known exact self-gravitating string solutions in a Minkowski 
background describe either infinite straight strings [3j HJ [5j [6] or infinite strings 
interacting with plane- fronted gravitational waves [71 [8], and the only general 
result available on the dynamics of string loops in strong-field gravity is Hawking's 
proof that a collapsing circular loop would form an event horizon after radiating 
away at most 29% of its original energy [9j. (In addition, the full 3-metric outside 
a circular loop at a moment of time symmetry has been constructed in [TU].) 

Given the complexity of the general strong-field problem, therefore, the most 
promising arena for the study of loop self-gravity would appear to be the weak- 
field limit. Fortunately, it can be shown that the self-gravity of a GUT string 
would almost everywhere be small enough to justify a weak-field treatment [TTj . 
although this treatment typically breaks down in the vicinity of certain common 
pathological features known as kinks and cusps. 

In [12] (henceforth referred to as Paper I) it was shown that a certain type 
of rigidly-rotating cosmic string loop - the Allen-Casper-Ottewill solution [13] 
evolves by self-similar shrinkage in the weak- field limit. This is the first explicit 
back-reaction calculation to have been published for a cosmic string loop, as all 
previous calculations have relied on numerical approximation [11] . However, the 
calculation relies on what might be called an adiabatic approximation, as the 
total flux of gravitational energy and momentum from the loop is integrated over 
a single oscillation period of the loop's unperturbed trajectory, and is shown to 
induce no changes other than a small decrease in the loop's energy and length, 
plus a comparable rotational phase shift. 

In this paper I extend the analysis of Paper I by relaxing the adiabatic as- 
sumption and constructing a self-similar evaporating solution of the linearized 
field equations which is identical to the Allen-Casper-Ottewill solution on spacelike 
sections and ultimately radiates away all its mass-energy. Although the trajectory 
of the string loop, and the perturbations it induces in the background metric, are 
only minimally different from the corresponding quantities in the adiabatic case, 
the evaporating solution represents an incremental advance in the understanding 
of loop self-gravity. In particular, the metric perturbations can be calculated at 
all points outside the string's worldsheet, and the perturbed (radiating) spacetime 
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can be matched to a flat spacetime which occupies the causal future of the final 
point of evaporation of the loop. The naive expectation that the loop evaporates 
completely to leave behind a vacuum remnant spacetime therefore appears to be 
consistent at this level of approximation. 

2. Preliminaries 

The mathematical formalism used to describe zero-thickness cosmic strings has 
been discussed in some detail in Paper I, and in this section I will offer only a 
brief summary. 

Throughout the paper, x a = [x°, x 1 , x 2 , x 3 } = [t, x,y, z] a are local coordinates 
on the four-dimensional background spacetime (M, g a b), and the metric tensor g a b 
has signature (+, — , — , — ), so that timelike vectors have positive norm. The world 
sheet T of the string is the two-dimensional surface it traces out as it moves, and 
is described parametrically by a set of equations of the form x a = X a (( A ), where 
the parameters ( A = ((°, C 1 ) are gauge coordinates. 

The intrinsic two-metric induced on T by a given choice of gauge coordinates 

is: 

lAB = 9abX a , A X b , B (2.1) 

where X a ,A is shorthand for dX a /d( A . The intrinsic metric ^ AB is assumed to 
be non-degenerate with signature (+, — ) almost everywhere on T. If 7 denotes 
I det(7^ B )| the Nambu-Goto action [UJ [15] has the form 

/ = -/*! 7 V2 rf 2 C (2.2) 

where fi is the (constant) mass per unit length of the string. 

The stress-energy tensor T ab induced by a given trajectory is found by varying 
the Nambu-Goto action with respect to the background metric g a b, according to 
the standard prescription 51 = — | J T ab 5g a bg 1 ^ 2 d 4: x, where g = \ det(g a ft)|. The 
equation for 51 inverts to give 

T a V) = "2<T 1/2 J (SC/5g ab ) 5 (4) (x c - X c )d 2 ( (2.3) 

where £ = —fi'j 1 ^ 2 is the Lagrangian density in (12.21) . and 

SC/Sgab = -\m ll2 l AB X\ A X b , B (2.4) 
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with 7 the inverse of 7^ B . 

In Paper I the gauge choice ^ AB = ki] ab was imposed from the outset, where 
Vab = diag(l,— 1) is the Minkowski 2-tensor. However, it turns out that this 
choice is not as convenient when the string is evaporating (and the trajectory 
is no longer time-periodic), and for this reason the more general versions of the 
string stress-energy and equation of motion will be retained in what follows. 

The extrinsic curvature tensor of the world sheet T has the general form 

K abc = p a m p bn V n p mc (2.5) 

where V„ is the derivative operator associated with g a b, and p ab = r y AB X a ,A X b ,B 
is the projection tensor corresponding to ^ab- In terms of the derivatives of the 
position function X a the curvature tensor K abc and its trace K° = g ao K abc can be 
written as 

R abc = qC^AC^BDxa^ ^ +p am {2 Q) 

and 

K c = q c d (l CD X d , CD +p mn T d mn ) (2.7) 

with q ab = g ab — p ab the orthogonal complement of p ab , and T^ c the Christoffel 
symbol associated with g ao . The equation of motion of the string is recovered by 
setting the functional derivative 5C/5X a equal to zero, and reads simply: 

K c = (2.8) 

In the absence of any sources of stress-energy other than the string itself, the 
system of equations for the metric and the string's trajectory is closed by imposing 
the Einstein equation G ab = -8nT ab , where G ab = R ab - \g ab R is the Einstein 
tensorQ Equations (12. 3p . (12.41) and (12.81) . together with the Einstein equation, con- 
stitute the strong-field back-reaction problem. As was mentioned above, the only 
known solutions to the strong-field back-reaction problem describe either infinite 
straight strings (with or without an envelope of cylindrical gravitational waves 
[31 HI [6]) or infinite strings interacting with plane-fronted gravitational waves 
(the so-called travelling- wave solutions [3 [S]). There are no known exact loop 

^^The convention adopted here for the Riemann tensor is that 

Rabcd ~2^9ac^bd ~^9bd^ac 9adibc 9boad ) 9ef i^ac^bd ^adJ^bc^' 

Also, geometrized units have been chosen, so that the gravitational constant G and the speed 
of light c are set to 1 . 
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solutions in a vacuum background, and an analysis of loop evolution is currently 
feasible only at the level of the weak-field approximation. 

In the weak-field formalism the field and trajectory variables are truncated at 
linear order in the mass per unit length \x of the string, which is assumed to be 
small compared to 1. This is certainly thought to be true of GUT strings, for 
which n ~ 1CT 6 . The underlying assumption is that the solution pair (g a b,X a ) to 
the strong-field back-reaction problem can be expanded as perturbative series of 
the form 

oo oo 

»o6 = 53A2 ) and X a = Y,^ k X a {k) (2.9) 

fe=0 k=0 
(k) 

where the functions g ab and Xfy are independent of /i. The weak-field back- 
reaction problem consists in setting equal to r] ab , the Minkowski metric tensor, 
and solving the Einstein equation and equation of motion (12.81) at linear order in 
fj, to obtain the metric perturbation g a l b and the trajectory functions XfL and 

In what follows, the Minkowski 4-tensor will be taken to have its rigid form 
i] ab =diag(l, —1, —1, —1). In solving the weak-field problem, a gauge must first 
be chosen for the metric perturbation g a l b \ and also for the functions XfL and X^. 

If h a b = fig^J the standard harmonic gauge conditions read h b a ,b= hh, a (where 
h = h b and indices are everywhere lowered and raised using rj ab and its inverse 
r] ab ). The Einstein equation G ab = —8iiT ab then reduces to Dh a b = — 16iiS a b, 
where Sab — T a b — \ rf lab^'c an d D = d% — V 2 is the flat-space d'Alembertian. 
With x a = [t,x\ a the retarded solution for h a b is: 

Mi,x) = -4 / d 3 x' (2.10) 

J |x-x'| 

where t' — t — |x — x'| is the retarded time at the source point x'. Since g 1 ^ 2 = 1 
to leading order in h ab , equations (I2.3P and (12.41) together give 



S ab (x' d )=fi J ^^{X^aX^b-^^X^aX^^ 

(2-11) 

with 7 1 / 2 and j AB evaluated using r] ab in place of g a b and X^ in place of X a . 
Furthermore, if g ab = r] ab then = to leading order, and at linear order 
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in h ab the equation of motion K c = reads: 

q^^X^cD -J AC 1 BD h ab X a , A X\ B X\ C D ) = -q cd l AB X\ A X\ B (h bd , a -\h aM ) 

(2.12) 

where q c d and 7 AC are evaluated using rj ab rather than g ab (but X a is retained in 
full). As was explained in Paper I, this equation is just the flat-space version of 
the Battye-Carter equation [TBI E] • 

If X a is now replaced by X^ + 5X a where 5X a = fJ-X^ then, with q c d and 
j AB evaluated using Xfo in place of X a , equation (I2.12p splits into two parts: 

t&f D XUcD=0 (2.13) 

and 

q c d7 CD 6X d , CD -2 W AB $X$ vA 5X a , B 7 OT Xf 0) , CD 

— —q cd 7 AB X( ) } AX(ty } B(hbd,a—^h a b,d) (2-14) 

at zeroth order and first order in /i respectively. (Here, round brackets on indices 
denote symmetrization, so that qaX?L, A = ^[q a X^, A +q^X^, A ].) 

The weak-field back-reaction problem, in its most general form, therefore re- 
duces to equations f[2T0]) . (12331) and f[2TT4D for h ab , X ( a 0) and 5X a . In order to 
specialize these equations to the problem in hand, it is necessary to first explain 
the choice that will be made for X^, and the corresponding choice of gauge 
coordinates (° and £ . This will be done in Sections 3 and 4 below. 

3. The Allen-Casper-Ottewill Loop 

The equation of motion at leading order (12.131) can be solved exactly by making a 
suitable choice of gauge coordinates. A common choice is the so-called standard 
gauge ( A = [t, <j] a , in which ^y AB is everywhere proportional to the Minkowski 
tensor rj AB =diag (1, —1). Since •y AB = r] ab X^, A X^,b the geometric significance 
of this gauge choice is that the tangent vectors X^ , T (which is timelike) and X? ^ , a 
(which is spacelike) are orthogonal and of equal magnitude. As was demonstrated 
in Paper I, the equation of motion (12.131) reduces in the standard gauge to the 
wave equation 
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If the subsidiary gauge alignment condition X? q Jt, a) = r is imposed then all 
loop solutions of (12.131) can be written in the form = [r, X(r, cr)] a , where 

X(r,a) = i[a(r + a) + b(r-a)] (3.2) 

for some choice of vector functions a and b. The gauge constraints r y TT + r y aa = 
and 7 CTT = together imply that |a'| 2 = |b'| 2 = 1, where a prime denotes 
differentiation with respect to the relevant argument. If the loop is in its centre- 
of-momentum frame then a and b are each periodic functions of their arguments 
with some parametric period L, and since X(r + L/2, a + L/2) = X(r, a) for any 
values of r and a when a and b are periodic, the fundamental oscillation period 
t p of the loop is L/2. 

In the centre-of-momentum frame, the total four-momentum of the string loop 
on any surface of constant t = r is: 

p a = H Xf 0) , r da = fi [1, |a'(r + a) + |b'(r - a)] da = 0] (3.3) 

Jo Jo 

and in particular the energy of the loop is E = fiL. The corresponding angular 
momentum of the loop is: 

J = lu I (axa' + bxb') da. (3.4) 



In the weak-field limit, the gravitational power P radiated by an oscillating 
compact source can be calculated from the standard expression for the power per 
unit solid angle in the wave zone (at distances large compared to the characteristic 
size L of the source): 

dP , ? 00 

^ = 7E m2 [ r tr^ fl f] (3.5) 

m=l 

where lj = 2^jt v is the circular frequency of the source, and 



f ab {m,n)=t; 1 [ T ab (t,x)e imu)(t - n -* ) dtd 3 x (3.6) 

Jo 

is the Fourier transform of its integrated stress-energy, with n the unit vector in 
the direction of the field point. 
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For a string loop in a flat background in the standard gauge, the stress-energy 
tensor (12.31) becomes 



T a \x c ) =fij (X ( a 0) , T X b (0) , T -Xfa„ X\ Q) „ ) 5^(x c - X ( c 0) ) dr da (3.7) 

and the power per unit solid angle reduces to 
dP 00 

^- = Siifi 2 m2 K A ■ A *)( B ■ B *) + \A-B*\ 2 -\A- B\ 2 }. (3.8) 



dn 

m=l 



where 



and 



A c {m,n) = {f e 27Tim[a+ ' n - a{a + )]/L [l,a'(a + )} c da + (3.9) 
Jo 

B c (m,n) = \ [ e 27Tim[(J — n - hi(7 - )]/L ll,b'((T_)] c dcT_. (3.10) 
Jo 

with the gauge coordinates a + and o"_ defined by a± = r±a. (Here and elsewhere 
the dot product is taken with respect to the flat-space metric r] ab .) 

The radiative efficiency 7 of a string loop is defined to be 7 = ^ 2 P, where 
P = J ^dVt is the total power of the loop. The radiative efficiency has been 
calculated analytically for a large class of simple loop trajectories by Allen, Casper 
and Ottewill [T5] , and numerically for the loop by-products of string network 
simulations by Allen and Shellard [18] and Allen and Casper [19]. The simulation 
studies found that the value of 7 is typically of order 65-70, and seems to be 
bounded below by about 40. This observation is consistent with the results of the 
analytic study, which found that 7 has a minimum value of about 39.0025, and 
occurs if the mode functions have the forms 

a(<r+) = {\a + \ - \L) z for - \L < a + < \L (3.11) 

and 



cr 



^[cos(2ttct_/L) x + sin(2vror_/L) y], (3.12) 
with a represented by its even periodic extension when a + lies outside [— |L, |^]J§ 



2 



Strictly speaking, the a mode function was represented in both [13] and Paper I in the form 



a(cr+) = 



<T+ i if < (7+ < \L 

(L-a + )z if \L < a + < L 
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The corresponding loop, which I will henceforth refer to as the Allen-Casper- 
Ottewill or ACO loop, is rigidly rotating about the z-axis with an angular speed 
uj = Atc/L. The evolution of the loop is illustrated in Figure 1, which shows the 
y-z projection of the loop at times r — e = 0, L/16, L/8 and 3L/16 (top row) and 
r — e = L/4, 5L/16, 3L/8 and 7L/16 (bottom row), where the time offset e is 
0.02L. The string has been artificially thickened for the sake of visibility, and the 
z-axis is also shown. The projections of the loop onto the x-y plane are circles 
of radius L/(4ir). The points at the extreme top and bottom of the loop, where 
the helical segments meet and the modal tangent vector a' is discontinuous, are 
technically known as kinks. They trace out circles in the planes z = L/8 and 
z = —L/8. All other points on the loop trace out identical circles, although with 
varying phase lags. The net angular momentum of the loop is J = ^[iL 2 z. 

I should stress that it is not currently known whether the ACO solution has 
the lowest radiative efficiency of all possible string loops. Allen, Casper and 
Ottewill [13] showed only that the ACO solution minimizes 7 over a special class 
of solutions with the mode function a given by (13. lip and the mode function 
b confined to the x-y plane but otherwise arbitrary. Nonetheless, in Paper I the 
first-order equation of motion (12.141) was solved for a leading-order trajectory 
of ACO form, and it was shown that the net perturbation 5X a induced in the 
string trajectory over a single oscillation period At = t p corresponds to no more 
than a change AL = —^j°fj,L ~ — 19.501/iL in the length and parametric period 
L of the loop, and a rotational phase shift with magnitude about 38.92/i acting 
to advance the overall pattern of the loop. In other words, as it radiates the ACO 
loop retains its shape and therefore its original radiative efficiency. If the ACO 
loop does indeed possess the lowest possible radiative efficiency, it follows that it 
would be the longest-lived loop in any ensemble of loops with the same energy. 

4. The Trajectory of the Evaporating ACO Loop 

Whether or not the ACO loop minimizes 7 , the fact that 7 remains constant 
as the loop evaporates is the key to a more extended description of the loop's 
evolution. Over the course of a double period At = 2t p = L the total energy 
E = fiL of the loop is reduced by an amount AE = —PL = — /z 2 7°L and so 

in which case the centroid of the loop lies at [x,y, z) = (0,0, \L). For the purposes of the 
analysis presented below it is more convenient to translate the loop down the z-axis by an 
amount hL so that the centroid lies at the origin. The resulting mode function is then given by 
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AE/At = — yU 2 7°, or equivalently AL/At = — fij°. Since the loop is rotating 
uniformly, this relation is also true instantaneously, so that dL/dt = — /i7° and 
therefore L(t) = Lq — /X7 (t — to), where Lq is the value of L at some initial time 
to- In particular, if the time coordinate is chosen so that to = then 

L(t)=L (l-t/t L ) (4.1) 

for t < ti, where tz = £o/(tV) * s the total future lifetime of the loop. 

If the adiabatic calculation performed in Paper I accurately represents the evo- 
lution of the ACO loop, it is expected that as it evaporates the loop retains the 
same spatial cross-section, apart from a uniform reduction in scale proportional 
to L(t), and a rigid rotation about the z-axis. What follows now is a heuristic 
argument leading to a plausible form for the trajectory of the evaporating loop. 
It should be cautioned that the argument below relies implicitly on an assump- 
tion that the background spacetime is flat, and consideration of the first-order 
equation of motion f!2.14p quickly shows that this assumption is overly optimistic. 
Nonetheless, it is also clear on physical grounds what form the relativistic correc- 
tions to the trajectory should have, and the first-order equation of motion will be 
solved to give these corrections explicitly in Section 6. 

The most general form of the position function X of the evaporating loop (with 
r henceforth replaced by the time coordinate t, and £ a general parameter) is 

X(U) = L(*)X(U) (4.2) 

where X is rigidly rotating, so that if t < t^ and t* < are two different times 
then 

cos — sin 
X(t, £) = sin0 cos0 
1 



X(t*,0 (4-3) 



with (f)(t*,t) a phase shift and £*(£,£*,£) a reparametrization of £. 

Here, if t* is fixed, X(t*, £*) is a reference function that at each value of t must 
trace out the same loop in M 3 as £ varies over its range. In particular, if t — > t + St 
then C(t + 5t,t*,£ + 5£) = £*(*, t*, f ) for some function <Sf (f). That is, 

dtC St + dt? 5^ = (4.4) 

to linear order in St and S£, and since this is a homogeneous linear first-order 
PDE with characteristic £ = a(t) for some function a, it follows that can be 
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any function of £ — ct(t). Hence, if the fixed reference time t* is suppressed, 



X(t,t)=L(t) 



cos0(t) — sin0(t) 
sin0(t) cos0(t) 
1 



X (£ - a(t)) 



(4.5) 



for some choice of functions 0, X and a. 

Furthermore, because the parametric period of a cosmic string loop scales with 
its length L, all timescales on the loop do likewise, and in particular the velocity 
field X t of the loop must at all times be rotating rigidly. Now, 



cos — sin 
sin cos 
1 



[L'Xo + L(0'MX o - a'X' Q )] 



(4.6) 



where 



M 



cos sin 
— sin cos 
1 



- sin — cos 
cos — sin 






" 


-1 


" 




1 




















(4.7) 



and the primes denote derivatives with respect to the relevant arguments. 

If the vector functions X , MX and Xq are linearly independent, the require- 
ment that X t be rigidly rotating constrains L', L<p' and La 1 to all be constant. 
The first condition, that V be constant, was derived independently at the begin- 
ning of this section, and is equivalent to requiring the radiative efficiency 7 to be 
constant. Since L(t) = L (l — tjti) for t < tL, the remaining constraints can be 
integrated immediately to give 



0(t) = -p\n(l-t/t L ) + 



and 



a{t) = -q ln(l - t/t L ) + a (4.8) 



where p, q, O and cto are constants. 

By suitably rotating the x-y coordinate axes and rezeroing the gauge coordi- 
nate £ it is always possible to set O and «o to zero. Furthermore, if a = p^/q 
then X is a function of a + pln(l — t/t L ) and so 



X(t,a) = L (l-t/t L ) 



cos0(t) — sin0(t) 
sin0(t) cos0(t) 
1 



X O (<7-0(O) 



(4.9) 



with 0(t) = -pln(l -t/t L ). 
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If this trajectory is to describe the evaporation of an ACO loop, the position 
function (14.91) and its time derivative should, to leading order in /i, match the 
position function X and velocity X 4 of the unperturbed ACO loop at time t = 0. 
The ACO loop is described by the mode functions (13.111) and (I3.12p . so its position 
function has the explicit form 

X{t,a) = ^L{cos[27r(t- ( r)/L]x + sin[27r(t-a)/L]y} + i(|t + or|-iL)z (4.10) 

for -\L <t + a <\L. 

Here, the constant scale factors Lq and L are taken to be identical. Since 
0(0) = 0, the position functions (14. 9 p and (I4.10p will agree at t = if 

X (a) = ^L[cos(27ra/L) x - S m(2na/L) f\ + l(\ a \-\L) z. (4.11) 

for — \L < a < \L. Furthermore, L'(t) = —L /tL is of order \i and can be ignored 
when calculating X 4 , while <p'(t) = p/ti at t = 0. Equating the time derivatives 
of (14.91) and (14.101) at t = then leads to the consistency conditions a = —4na/L 
and p = 4ntL/L = 4ti / (j jj) . After substituting these conditions and equation 
(14.111) into (14.91) . the position function of the evaporating loop becomes: 

X{t,a) = ^(l-t/t^lcosti^ + ^x + sinti^ + ^y + itl^^-al-Trlz} 

(4.12) 

for -2tt < (f>(t) - a < 2vr, where (j){t) = -{±irt L /L ) ln(l - t/t L ). 

Equation (14.121) suggests that a convenient choice of gauge coordinates (C°, C 1 ) = 
(u, v) for the trajectory of the evaporating ACO loop is 

u = \{(j){t) - a) and v = ±(<f>(t) + d). (4.13) 

Then 1 - t/t L = e ~ K ^( u + v ) ) where K[i = (Airt L / Lq)^ 1 or equivalently k = 7°/(47r). 
This in turn means that the time coordinate t on the loop can be written as t = 
t L [l — e ~ K ^( u + v )^ anc l the world sheet T of the evaporating loop can be represented 
by the 4-function 

X>, v) = [t L , Of + ±L e-****) [-1 r(u, v)] a (4.14) 

where 

r(u, v) = cos(f ) x + sin(t>) y + (\u\ — \n) i (4-15) 
for -7r < u < 7r, and t L = t^; l o- 
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Strictly speaking, the vector function r appearing in the 4-function (I4.14p is 
the 27r-periodic extension (in u) of the function defined in (14.151) . The gauge 
coordinate pairs (u, v) and {u + 2nk, v — 2nk) therefore correspond to the same 
spacetime point for any integer k. However, for the sake of simplicity the range 
of the parameter u will henceforth be restricted to (— 7r, tt]. The final evaporation 
point x a = [tL, 0} a of the loop then corresponds to the limit v — > oo. Note that in 
this limit the scale factor L(t) = L e~ K ^ u+v > goes to zero, while the angular speed 
dv/dt of the loop diverges and r(u,v) is undefined (although the loop's 4- velocity 
X t a , being scale-free, remains subluminal everywhere). 

If \nfi(u + v)\ <C 1 the loop's position function (I4.14p can be expanded in 
powers of \i to give 

Xfo{u,v) = ±L [u + v,r(u,v)] a (4.16) 

and 

5X a {u, v) = -±L k^{u + v) 2 , {u + v) r{u, v)} a . (4.17) 

It is these functions that need to satisfy the zeroth- and first-order equations of 
motion (I2.13P and (I2.14p . Note that the leading-order position function (14.161) is 
effectively generated by an expansion of X a about the curve u + v = on the 
loop's world sheet, or equivalently about the t = spacelike cross-section. 

The position function X a can be expanded about other cross-sections of con- 
stant t < tL by writing v — + v where v* is any fixed real number. Then 
equation (14.161) becomes 

X a (u,v) = [t L ,0} a + ±L*e- K ^ u+ ^[-±r(u,v*+v)} a (4.18) 

with = Lq e~ KfJ - v * , and if \k^l{u + v)\ <C 1 the zeroth-order position function is 

X? 0) (u,v) = [t^0} a + ±L4u + v,v(u,v* + v)} a (4.19) 

where = ^^L (l — e~ K/iV *). The analogue of equation (14.171) for 5X a is found 
by replacing L with L*, u + v with u + v, and r(w, v) with r(u, i>* + v). Clearly, 
any two expansions of this type will be related by a time translation, a constant 
dilation, and a rotation in the x-y plane, and if any of the expansions satisfies the 
equations of motion (12.131) and (I2.14p then they all will. 

The equations of motion will be considered more closely in Section 6. It will 
come as no surprise that the zeroth-order position function ( 14.161) does satisfy the 
zeroth-order equation of motion (12.131) . which reads simply X^, uv = 0. Before 
considering the first-order equation of motion (12.141) . however, it is necessary to 
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calculate the metric perturbations h a b, which depend only on X? y This will be 
done in Section 5. 

After a long calculation described in detail in Section 6, it turns out that the 
first-order position function 5X a given by (I4.17P does not satisfy the first-order 
equation of motion (12.141) . The reason for this is not difficult to appreciate. As 
was mentioned at the beginning of this section, the calculation leading to the 
evaporating loop trajectory (14.141) was predicated on the assumption that the 
background metric is flat. It is evident from (I2.14p . however, that the first-order 
equation of motion involves the first-order metric perturbation h^. Once the 
background metric g a b ceases to be flat, the spacetime coordinates [t,x,y,z] no 
longer measure proper times or proper lengths. In particular, the self-similar 
shrinkage driven by the scale factor L e~ K ^ u+v ^ may proceed at different rates in 
different directions at various points on the world sheet T. 

This possibility is easily accounted for by adding scale corrections, of order /j,, 
to the components of the evaporating loop's position function (14.141) . The form of 
these scale corrections is of course constrained by the equation of motion (I2.14p . 
but even so there remains considerable freedom in the choice of the correction 
functions, corresponding not only to possible transformations (at order /i) of the 
gauge coordinates (u,v) and the spacetime coordinates [t,x,y,z], but also to the 
propagation of free vibrations (with amplitude fi) around the loop. 

In what follows, the freedom in the choice of gauge coordinates will be partially 
fixed by taking the position vector to have the form 

X a (u,v) = [t L ,0} a + ±-L e- K ^ u+ ^[-^ 2F ^ v \f(u,v)] a , (4.20) 

where now 

r( u , v) = e^'^cos^) x + sm(v) y] + [(\u\ - ~tt) + fxA(u, v)} z (4.21) 

for — 7r < u < 7i. Here, the correction functions F, G and A are assumed to be 
continuous and 2-7r-periodic in both u and v , so that (u, v) and (u + 2tt, v — 2tt) 
again correspond to the same spacetime point, and all three functions remain 
small compared to as u + v — » oo. The assumption that the trajectory is 
dilated by the same factor e^ G in the x and y directions effectively removes any 
freedom to redefine the gauge coordinate v% 

3 To see this, suppose that the term G(u, v){cos(v) x + sin(i>)y} in equation (|4.23[) has the 
more general form G{u, v) cos(u) x + H (u, v) sin(v) y- Then if v is replaced by a new coordinate 
v' defined by v = v' + [J.V(u, v 1 ) for some function V, the correction term G cos(-u) ii + H sin(u) y 
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The presence of free vibrations which break the self-similarity of the world- 
sheet can be removed by requiring the intrinsic 2-metric ^ AB = g a t,X a ,AX b ,B to 
have the strictly self-similar form 

7ab(«, «) = (^L ) 2 e~ 2 ^ + ^ AB (u) (4.22) 

to linear order in /i, with ^ AB a 2-metric to be determined. This constraint will 
be considered in more detail in Section 6.2. 

Given ( 14.20(1 . the leading-order position function XZ* near t = retains the 
form (14.161) . while the first-order position function becomes 

5X a {u,v) = -±L K^{u + v) 2 ,{u + v)r{u,v)] a 

+ ^L fi[~F(u, v), G(u, v){cos(v) x + sia(v) y}+A(u, v) z] a . 

(4.23) 

The actual form of the correction functions F, G and A will be determined once 
the first-order equation of motion ( 12. 14ft has been solved, and additional gauge 
conditions have been imposed. 

Two remarks about the final form (14.201) of the loop trajectory are in order 
here. The first is that although (14.201) was developed to describe an evaporating 
ACO loop over the time interval < t < ti (or equivalently < u + v < oo), it 
can without difficulty be assumed to extend over the entire interval — oo < t < th 
(or — oo < u + v < oo). The trajectory then describes a loop that has been 
evaporating "eternally", with its energy and length scaling (to leading order in 
jj) as 1 — tjth- This is obviously just a mathematical idealization, but it does 
side-step the need to construct initial data for the metric perturbations h a b at 
t = 0. 

The second point is that the value of the constant k appearing in the corrected 
trajectory (I4.20p is initially undetermined. Although it is to be expected in view 
of the results of Paper I that k ~ 39.0025/(47r) w 3.1037, no particular value of 
k is assumed in the analysis that follows. It is only after the first-order equation 
of motion is solved that it becomes evident that the periodicity constraint on G 
forces k to take on its expected value. 

A schematic representation of the evaporation of the loop is shown in Figure 
2. The thickened line corresponds to the outer envelope of the loop, whose radius 

becomes (G cos v' — V sin v') x + (H sin v' + V cos v')y to leading order in ^. The unique gauge 
choice V — (G — H ) sin v' cos v' then reduces this term to the isotropic form G'{cos(V)x + 
sin(w') y}, with G' = Gcos 2 v' + H sin 2 v'. 



15 



shrinks to zero at the final evaporation point x a = [tL,0] a . Also shown are the 
future light cones Fq and Fl of the origin and the evaporation point respectively. 
Note that, because of the first-order corrections to the background metric and 
the loop trajectory X a , the envelope of the loop and the light cones depicted in 
Figure 2 should deviate from straight lines by terms of order \x. 

Now that the choice of the evaporating loop trajectory f )4.20p has been ex- 
plained in some depth, the remainder of the paper is a straightforward applica- 
tion of the solution methods discussed in Section 2. The metric perturbations 
h ab generated by X? s are calculated in Section 5, the equation of motion for the 
perturbation functions SX a is constructed and solved in Section 6, and the metric 
perturbations on the future light cone Fl of the final evaporation point are shown 
to match onto an empty, flat spacetime in Section 7. 

5. The Metric Perturbations 

5.1. Rotating self-similarity of the metric perturbations 

The metric perturbations h ab at any field point x a = [t,x] a away from the world 
sheet of the evaporating loop can be calculated from equation (I2.10p . Because 
the analysis of Section 7 depends on an understanding of the behavior of h a b 
arbitrarily close to the light cone Fl, it will be necessary to consider source points 
on the evaporating loop close to the final evaporation point (where u + v —>■ oo) 
and so it cannot be assumed that \kh(u + v)\ <C 1. The scale factor Lq e,- K ^ u+v ) 
will therefore be retained in full in this section, although the correction functions 
F, G and A - which are bounded and appear only at order \i in (j4.20f) - will be 
ignored. In other words, the trajectory will in this section be assumed to have its 
flat-space form (14.141) . 

The first step in the calculation is to write down the leading-order expressions 
for the intrinsic 2-metric 7 AB and the source function S a b- Since 

X a , u = ±L e -M«-H0 [i, sgn ( u ) g]« (5.1) 

and 

X\ v = ±L e -«A«(«+*) [1, - sin(^) x + cos(v) y} a (5.2) 

to leading order in Li, it follows that the intrinsic 2-metric 7 AB = X,a -X, b on the 
world sheet has components 

luu = l m = and luv = (±L ) 2 e- 2 ^ u+ ^ (5.3) 
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to the same order. Hence, 7 1 / 2 = (-^L ) 2 e 2k K u + v ) ( while the inverse 2-metric 

Ml 



7 has components 



The source function defined in equation (12. lip can be written as 



S ab (x' d )=fi / y ab (u,v)5 (4 \x' d -X d )dudv 



(5.5) 



where 

V ab = 7 1 /V B (^a,A -X, B ), (5.6) 

evaluated again at leading order in \i. Note here that ^ AB X ai A X b ,B is just the 
world-sheet projection tensor p a &, which in matrix form is 



Pab 



2 sint> — cosv —s 

sinv — ssinv 

-cosv scosf 

— s — ssint> scosf 



(5.7) 



J ab 



where s = sgn(u). The function ^! ab is therefore equal to — 7 1 ' /2 q , a b, where q ab = 
Vab~Pab was introduced in Section 2 as the orthogonal complement of p ab . Hence, 
y ab = (±L ) 2 e- 2 ^ u+v ^ ab with 



ab{S,V) 



1 

sinv 

- COSf 

— s 



smv —cosv —s 
1 — ssinv 

Is cos v 

-ssint> scosf 1 



(5i 



ab 



Given (15.51) . the integral expression (I2.10p for h ab reads 



h ab (t,x) = -4/i J ^ - — J 



y ab {u,v)5 (4 \x' d -X d )dudv (5.9) 



where the source point is x' d = [t',x'] d , with t' = t — |x — x'| the retarded time^. 
The distributional factor 5 (4) (x /d - X d ) in ([53]) can be integrated out by first 



4 Note that, at a physical level, the equation t' = t — |x — x'| imposes the requirement that 
the source point [t',x'] lie on the backwards light cone of the field point [t,x] to leading order 
in \i. 
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transforming from y a = [u, x'] a to z a = [t^x'] — X a (u, v), with v, t and x held 
constant. The Jacobian of this transformation is: 

\dy a /dz b \ = \dz h /dy a \~ l = \X°, u -n-X, u \-\ (5.10) 

where n = (x — x')/|x — x'| is the unit vector from the source point x' to the field 
point x. 

On integrating over z a , t' is everywhere replaced by X° = ti[l — e~ K ^ u+v ">] 
and x' by X = -^Lq e" K ^ u+v 'r(u, v), while u becomes an implicit function of v, 
t and x through the equation X°(u,v) = t — |x — X(/u, v)\. The integral for h a b 
therefore reads 

M*,x) = -4/i J ||x-X|X°, u -(x-X)-X, u |- 1 ^ afc ( M ,t;)^. (5.11) 

where, as before, X°, u = ^L e -^i u + v ) an d X, u = -^L e~ K ^ u+v ^sz to leading 
order in \i. 

Evaluation of the integral (15.111) is simplified by the fact that h ab satisfies a 
very elegant self-similarity relation. If [t, x] is any field point with t — ti < |x| 
(where t — t L = |x| is just the equation, at leading order, of the future light cone 
Fl of the final evaporation point), first define 

^{t, x) = -(nfi)- 1 ln[(|x| - t + t L )/t L \ (5.12) 

then consider a second field point [i, x] with 

t = e K ^'|x| and x = e^i2(--0)x (5.13) 



where 

cos 9 — sin 9 

R{9)= sin# cosfl 

1 



(5.14) 



is the standard rotation matrix about the z-axis. Note that [t, x] lies on the future 
light cone F of the origin, as t = |x|. 

Now, the equation t — X°(u, v) — |x — X(m, v )| which specifies u as a function 
of v in the integral (15.111) - or equivalently on the backwards light cone of the 
original field point [t, x] - reads explicitly 

t-t L [l- e -^ iu+v) ] = |x - ±L e - K ^ u+v) t(u, v) I . (5. 15) 
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Multiplying this equation by e K/ ^(*' x ) gives: 

(t - t L )e K ^ + t L e~ K ^ u+v -^ = \e K ^ x - ±L e -«f*(«+»-i« r ^ ^ | 

where (t — t L )e K ^ = e K ^|x| — t L and 

r(u,v) = cos(f ) x + sin(t>) y + (\u\ — \ti) z 

= R(ip)[cos(v — ifj) x + sin(t> — ip) y + (|w| — |7r) z] 
= R(ip) r(u,v — ip). 



(5.16) 



(5.17) 



Since |v| = |-R(— ip)v\ for any 3-vector v, and R(—tp)R(tfj) = I, equation (15.161) 
can be rewritten as 



t-t L [l-e 



-Kfl(u+V — llj)'l 



lL e- K ^»r(«,.-^)|, 



(5.18) 



or equivalently t — X°(u, v — ip) = |x — X(w, v — ip)\. 

Hence, if u = U (v) is the equation of the curve of intersection of the world sheet 
T with the backwards light cone of a given field point [£,x], then the equation 
for u on the backwards light cone of the corresponding image point [t, x] on F is 
u = U(v — ij)(t, x)). Similarly, the integrand in equation (15. lip for h a b, which is 
equivalent to 



^L n e~ Kfl{u+v) \\x-X\ 



■In 



u,v) 



S [X - X(U,V)] • Z I ^abis, V) 



(5.19) 



can be re-expressed as^l 

f L e -^C«+«-^) | | e K ^ J R(-^)x - X{u, v 



„• • - - S [e^H)x-X( U> « 

x^ afe (s,f) 

e -*f («+*-^) | | x - X(u, v-ip)\ - s [x - X(u, v - V)] • z I" 1 

d( 



xR a c (^ cd (s,v-m b d W 



where 



R:^) = (i®Rm: 



10 

cos ip — sin ip 

smip cos ip 

1 



(5.21) 



5 Note here that 



6 k^ j x _ X ( U; „)] . 2 _ [ e K ^i?(-</;)x - X(u, u - V)] • z 



for the simple reason that the z-component of x is unaffected by the rotation matrix R, while 
the z-component of X is ±L e - K ^ u+v) (\u\ - 



(5.20) 
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Comparison of equations (15.191) and (15.201) . together with the fact that if u = 
U{v) on the backwards light cone of [t, x] then u — U(v — if) on the backwards 
light cone of [t, x] , indicates that 

h^t, x) = R^) h cd (i, x) R b d {iP), (5.22) 

where ip(t,x) is given by equation (15.121) . The geometric meaning of this relation 
is illustrated in Figure 3. To calculate h ab at a field point [t, x] in the past of F L , 
simply dilate the point by a factor e KfJj ^ through [t^, 0] and rotate it by an angle 
—ip in the x-y plane to map it to the image point [i, x] on Fq. Then calculate h a b 
at [i, x] and rotate each of the components of the resulting tensor by an angle ifi 
in the x-y plane to give h a b at [t, x]. 

From a computational viewpoint, equation (15.221) allows us to calculate h a b 
at any point away from the world sheet of the loop from a knowledge of h a b on 
the null surface Fq alone. Furthermore, it points to the existence of a symmetry 
in the background metric induced by the self-similarity of the loop's trajectory. 
This symmetry can be better appreciated by writing the components ( I5.13P of the 
image point explicitly in terms of t and x: 

t = t £ (l + |w|) -1 |w| and x = t L (l + |w|)~ 1 i?(-^)w (5.23) 

where 

w = x/(t L -t). (5.24) 

If it were not for the presence of the rotation operators in fl 5 . 2 2 1) and (I5.23p . 
which depend on if) rather than w, the metric r\ ah + h a b would be a function of w 
alone and so would be strictly self-similar. As it is, the rotation operators act to 
enforce a slightly more complicated symmetry in the metric which could perhaps 
be termed "rotating self-similarity" . 

5.2. Explicit form of the metric perturbations on Fq 

Suppose now that [|x|,x] is a general point on the future light cone F of the 
origin. The equation (I5.15P for function of v on the backwards light cone 

of [|x|, x] reads: 

|x| - t L [l - e~^ iu+ ^} = |x - ±L e~ K ^ u+v) r(u, v) \ . (5.25) 

Application of the triangle inequality to the norm on the right-hand side of this 
equation then gives: 

t L \l-e-^ {u+v) \ < j-L G e- K ^ u+v) [l + (|«| - |tt) 2 ] 1/2 , (5.26) 
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as \r(u, v)\ = [1 + (\u\ — l^) 2 ] 1 / 2 . Since ti = ^t^Lq and the maximum value of 

[1 + (\u\- \ti) 2 } 1/2 for u G (-7T, tt] is (1 + ^vr 2 ) 1 / 2 « 1. 8621, this inequality can be 
re-expressed as 



| e «^(«4*) _ !| < (! + l n 2 ) 1/2 Kfi. (5.27) 

It is therefore evident that \nfi(u + i>)| will remain small at all points on 
the backwards light cone of [|x|,x] provided that Kfi itself is sufficiently small. 
For example, \nfi(u + v)\ will be everywhere no greater than 10~ 2 if ln[l — (1 + 
k/j] > — 10 -2 , and this in turn will be true if Kfi < 5.34 x 10~ 3 . Given 



that we expect that k ~ 3.1037 and \x will be of order 1CT 6 or smaller for a GUT 
string, the assumption that \k/j,(u + v)\ remains small is not an unreasonable one. 

Expanding e ~ K f 1 ( u + v ) ' m powers of fi in (15.251) and retaining only the leading- 
order terms gives the following equation for u as a function of v. 

T — (u + v) = [(X — cost;) 2 + (Y- sinv) 2 + {Z - (\u\ - ^vr)} 2 ] 1 / 2 , (5.28) 

where 

(X,Y,Z) = 4vrL 1 (x,y,^) and T = (X 2 + Y 2 + Z 2 ) 1 ' 2 . (5.29) 

Equation (I5.28j) can easily be solved to give an explicit formula for u, but it turns 
out that in calculating h a b it is more useful to regard u as the independent variable. 

Similarly, for small values of \k/j,(u + v)\ the integrand ( 15. 19ft in the formula 
( 15.111) for hab at [|x|,x] becomes: 

{[(X-cost;) 2 + (r-sin W ) 2 + {Z-(| M |-i7r)} 2 ] 1 / 2 - S [Z-(|n|-i7r)]}- 1 ^( S ,t;). 

(5.30) 

Note that the denominator in this expression is manifestly positive or zero. In view 
of equation (15.281) and the fact that s\u\ = u, the denominator can be expressed 
in the simpler form 

[(X-cosv) 2 +(Y-smv) 2 +{Z-(H-^)} 2 ] 1/2 -s [Z-{\u\-^)] = T-s(Z+\it) -z 

(5.31) 

and (15. lip reduces to 



MI*|,x) = -4// J[T - s{Z + \ii)-v]- 1 ^> ab {s,v)dv. (5.32) 

All that remains now is to determine the limits of integration. The entire curve 
of intersection V of the world sheet with the backwards light cone of [|x|,x] will 
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be traced out as u varies from — 7r to ir. The integral (I5.32p therefore breaks into 
two parts: one with s — — 1 and u G [— n, 0], and one with s = +1 and u G [0, 7r]. 
The corresponding values of v, which I will denote as Vk = v\ u= k-n for k = —1,0 
and 1, are roots of the light-cone equation (15.281) : 



T — (kit + Vk) = [(X-cosV k ) 2 + {Y - sinV k ) 2 + {Z - (\k\ - |)tt} 2 ] 1/2 . (5.33) 



Although it is not possible to solve for Vk explicitly except in certain special cases, 
it is clear by inspection that V-i = V\ + 2tt. Furthermore, implicit differentiation 
of (I5.28P with [T, X, Y, Z) constant indicates that 



Some of the mathematical consequences of equations (I5.33P and (15.341) are 
examined in Appendix A. In particular, it is shown that the values of Vq and V\ 
are bounded on Fq, and that v is a monotonic decreasing function of u. The last 
property follows partly from the fact that the numerator in (15.341) is positive or 
zero from (15.311) . while the denominator is positive or zero as a consequence of 
fl5T28l) and the identity 

(X — cosv ) 2 + (Y — sin^) 2 = (X sin v — Y cos v ) 2 + (X cosv + Y sin v — l) 2 . (5.35) 

However, cases where either the numerator or denominator is zero for some 
value of u need to be considered separately, and it is shown in Appendices A. 2 
and A. 3 that the corresponding spacetime points have interesting physical inter- 
pretations. Perhaps the most important result is that the spacetime derivatives of 
hgh are undefined at all points with the property that the denominator in (15.341) 
is zero at one of the extreme values u = or ±7r. 

Given that VLi = Vi + 2tt > Vq > Vi, the equation for the metric perturbations 
becomes: 



dv 
du 



T-s(Z + \n)-v 



(5.34) 



T — {u + v) + X sin v — Y cos v 



M|x|,x) = -4/i / [T + (Z + ±n)-v}- 1 y ab (-l,v)dv 




JV 
rV 



-4// / [T — (Z + |tt) - ^]- 1 ^ a6 (+l, v) dv. 



JVi 



(5.36) 



After referring back to the definition (I5.8P of \l/ a b, it is readily seen that 





/+ + /_ S + + S_ -C + -C_ 
S + + S_ /+ + /_ S + -S_ 




(5.37) 



I+-I- S+-S- -C+ + C- + F 
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where 

I+ = H(x + -Vi-2k)/(x + -V )], 



- s + - 









and 



- cos x+ sin 
sin x+ cos x- 



— cosx- sm%_ 
sin x„ cos x_ 



I_ = ln[(x_-V )/(x_-V 1 )], 

Si( X+ -Vx~ 2n) - Si( X+ - Vb) 
Ci( X+ - VI - 2tt) - Ci( X+ - Vb) 



(5.38) 
(5.39) 



Si( X _ - Vb) - Si( X _ - VI) 
Ci( X _ - Vb) - Ci( X _ - VI) 



(5.40) 
f w~ x cos if tito 

J X 



with x ± = T±(Z+|), and Si(x) = sin w and Ci(x) = 

the usual sine and cosine integrals. 

One case in which it is possible to solve explicitly for Vb and Vi, and so to 
write explicit formulas for the metric perturbations h a b, occurs when the image 
point [|x|,x] lies on the z-axis. Then X = Y = and T = \Z\, and equation 
( 15.331) reduces to 



V k = \Z\ -kn-[l + {Z 



(5.41) 



In particular, 



V = \Z\ -[1 + {Z+ I) 2 } 112 and V x = \Z\ - tt - [1 + {Z - f f} 1 ' 2 (5.42) 

while X ± = \Z\ ±{Z + \). 
Hence, 

( X+ - VI - 2n)/( X+ - Vb) 

= (x--v )/(x_-v 1 ) 

= {Z-\ + [l + (Z-f) 2 ] 1/2 }{-^-f + [l + (^ + f) 2 ] 1/2 } (5-43) 



and 7_ = J+. This means that ht z = 0, and that 

2 L V 2- 



8^1n({Z - | + [1 + (Z - \) 2 \ ll2 }{-Z - § + [1 + (Z+ I) 2 } 1 ' 2 }) 



(5.44) 



increases monotonically from approximately — 19.734/i (at Z = 0) to (as |Z| 
oo). In fact, for large values of \Z\, 



h tt 



h 



h, 



\kh/\z\ + o{\z\ 



-2\ 



(5.45) 
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The variation of hu/iAfi) with Z is plotted in Figure 4a. 

Explicit formulas for the remaining non-zero metric perturbations, h tx = 
4/i(5+ + S-), h ty = -4 M (C+ + C_), hzx = 4/x(S+ - SL) and = -4/i(C + - C_) 
can be found by substituting the expressions (15.421) into (I5.39j) and (15.401) . How- 
ever, there is little to be gained from writing these formulas out in full. It is 
easily seen that if Z is replaced by — Z then S+ «-> — S- and C + <-> — C_, and so 
htx ~^ —htx, h ty — > —hty, h zx — > and /i zy — * h zy . This set of symmetries is to 
be expected on physical grounds, as the ACO loop is invariant under a complete 
spatial inversion (x — > — x, y — > — y and 2 — > —z), and therefore in this special 
case the space-space components of h a b should be even functions of Z, and the 
time-space components odd functions of Z. 

The variations of h tx /(4:fi), h ty /(4:fi), h zx /{Aii) and /i 2J/ /(4/i) are plotted in 
Figures 4b to 4e respectively. At Z — 0, the perturbations h tx and h ty are of 
course both zero, while h zx w —5. 308/i and /i 2J/ w 12. 395/i. All the perturbations 
tend to zero as \Z\ — >• oo, with 

^ = 47r /U /Z + 0(|Z|- 3 ), (5.46) 

= -4tt/V(Z|Z|) + 0(|^r 4 ), (5.47) 

h zx = -Anfi/\Z\ +0(\Z\- 3 ) (5.48) 

and 

h zy = f7r/i/Z 4 + 0(|Z|- 6 ). (5.49) 

6. Solving the String Equations of Motion 
6.1. Near- field form of the metric perturbations 

For the purposes of setting out the first-order equation of motion (12.141) . it is neces- 
sary to generate expressions for the metric perturbations h a b in the neighborhood 
of the world sheet of the evaporating loop. As was mentioned in Section 4, when 
doing this it is sufficient to restrict attention to points near the spacelike surface 
t — only. So consider a field point x a = [t, x] a of the form X a (u, v) + 5x a , where 
X a {u, v) is a fixed point on the flat-space trajectory (14. 14j) with \k/i(u + v)\ 1, 
and the components of the displacement 6x a are all of the order of some small 
parameter. 

A first point to note is that if X a (u, v) is any source point on the loop's 
trajectory on the backwards light cone of X a (u,v) then \k[i{u + v)\ will also be 
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small compared to 1 for reasonable values of Kfi. This follows from equation 
( 15.151) . which in the current situation reads: 

t L [1 - e -»"( u +»)] -t L [1 - e -^{u+v)^ = | ^ Lq e -^(u+v) r ^ ^ _ e -« M (a+e) ^ | _ 

(6.1) 

An application of the triangle inequality gives 
f£ | e -*M(fi+e)_ e -«M(«+«)| < j_ Lo e -M«4«) | r ( UjV )| + J_^ oe -«M(o+«) | r (w,u)|, (6.2) 

or equivalently (in parallel with the derivation of (I5.27P ). 



e KH(u+v)-K/J,(u+v) I < ^ _|_ e /t/i(u+u)-«/x(u+«)^-^ _|_ ly]- 2 ) 1 / 2 /-^. (6.3) 



Hence, 



+ f ) — /t/i(-u + v)\ < In 



1 + (1 + ^TT 2 ) 1 / 2 ^ 

1 - (1 + V) 1 ^ 



(6.4) 



and so, for example, + u)| will be no greater than 10 everywhere if ft// < 

2.69 x 1(T 3 . 

Since \k^{u + v)\ is guaranteed to be small, the equations (I5.33j) for and 
f l5.37l) - f)5.40p for /i a6 continue to apply in this case, with 

T = u + v + 4vrLo 1 5t and (X, Y, Z) T = r(u, v) + 4ttL x 5x. (6.5) 

where as before r(w, v) = cos(f) x+sin(t>) y + (|w| — z. In particular, if 5x a = 
the equation for becomes: 

u + v - {kit + V h ) = [2 - 2 cos(t; - Vfc) + (H - |fc|vr) 2 ] 1/2 , (6.6) 

and, given that u is restricted to the range (— 7r,7r], the trivial solution Vfc = t> 
holds in two cases: (i) if u > and = 0; and (ii) if u < and k = —1. 

In the analysis that follows it will be assumed that u > 0. Some remarks 
about the behavior of h a b near points X a (u,v) with it < will be made towards 
the end of this section. When u > the z-component of the field point satisfies 
Z + = u + 4ttLq Sz, and so 

X± = v+ (u±u) +4vrL 1 (5t ±fc). (6.7) 

Furthermore, if 5x a = the limits of integration in (15.361) are given by Vq = v 
and Vi — v — W(u), where W is the unique root of the equation 

u - n + W =[2-2cosW + (u-n) 2 ] 1/2 . (6.8) 
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Note that because V > V\ it follows that W > 0, while it is clear from (16.81) that 
W = only at the kink point u = it. For future reference, it is useful to rearrange 
(16. 8p in the form 

u - ir = W~\l - cos W) - \W. (6.9) 

(Incidentally, one way to verify that (16.81) has a unique solution W is to note that 
the function on the right-hand side of (I6.9P is strictly decreasing.) 
Implicit differentiation of the defining equation (15.331) for V k gives 

d T V k = H k 1 [T-(kn + V k )} 

and 

(d x V k , dyV k ,d z V k ) = H k \cosV k - X,sinV k - Y, (\k\ - |)tt - Z), (6.10) 

with 

H k = T— (kir + V k ) +XsmV k -YcoaV k . (6.11) 

Note that if 5x a = then H = u and Hi = u — tx + W — sin W. 
Hence, to linear order in 5x a , 

V = v + AnL^iSt - 8z) and Vi=v-W + A-kLq l 8V (6. 12) 

where 

5V = H{ 1 [(u-n+W)5t+{cos(v-W)-cosv}5x+{sm(v-W)-smv}5y+(n-u)5z}. 

(6-13) 

Unfortunately, the linear expansion (16.121) for Vq is not sharp enough to resolve 
the singular behavior of h a b near the world sheet, as \- = v + AttLq x (5t — 8z) and 
so X- ~ = at this order. 

At quadratic order in 5x a it turns out that 

V = v + AitL Q \5t - 8z) - \u-\AirL l ) 2 Q + , (6.14) 

where 

Q + = (8x) 2 + (8y) 2 + (8t - 8z)(8t - 8z + 25xsinw -25ycosv) 

= ^ ab (+l,v)5x a 5x b . (6.15) 

Since ^! a b is proportional to the orthogonal complement q a b of the projection tensor 
p a b, the quadratic form Q + is zero if and only if 5x a lies in the tangent plane to 
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the world sheet at X a (u,v). It should therefore be understood that in evaluating 
h ab and its derivatives the limit Sx a — > is taken only from directions outside the 
world sheet. 

The algebraic computations whose results appear throughout the rest of this 
section are long and very mechanical, and their details have been relegated to 
Appendix B. Neglecting contributions that tend to zero in the limit as 5x a — > 0, 
the metric perturbations have the near-field form: 



{^Kb = {ln[(4 7 rL 1 ) 2 Q + ] -\n(2u)}V ab (+l,v) + ( 1e - CiW) $ ab (+l,v) 

+ (Si W) n ab (+l, v) + K{u) $ a6 (-l, 2u + v) - E(u) n ab (-l, 2u + v) 
-(In W) n ab (+l) + [In (W + 2u - 2tc) - ln(2«)] n a6 (-l), (6.16) 

where 7 E xs 0.5772 is Euler's constant, 



®ab(s,v) 



sin v — cos v 

sint> — ssin-y 

— cosf scosf 

— ssinf scosf 

Q ab (s,v) = d v $ ab (s,v) 



- ab 



and 



ru(s) 



1 -s 

10 

10 

-s 1 



(6.17) 
(6.18) 

(6.19) 



ab 



Note in particular that ty ab (s,v) = § ab (s,v) + U ab (s). Also, the functions K and 
E are shorthand for 

K{u) = Gi{W + 2u-2ti)-Ci{2u) and E(u) = §i{W + 2u-2n) -Si (2m) . 

(6.20) 

Similarly, the leading-order contributions to the spacetime derivatives of h ab 

are: 



{A^)-\j-L Q )d t h ab 



(±L ) {d t A )#„ & (+l, v) + {St - Sz) (d t A ) n ab (+l, v) 
+ [A + 7e - Ci W - ln(2u)] il ab (+l, v) - (SiW) $ a6 (+l, v) 
fi„ b (-l, 2m + u) + E(u) $ a6 (-l, 2m + u) 



sin 



u-jr+W-siaW W 
1 



[± v-W) 



W+2u-2w 



t> ab (-l,v-W)} : 



(6.21) 
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{Aii)- x {±L )d x h ab = (±L )(d x A)V ab (+l,v) + (5t-5z)(d x A)n ab (+l,v) 

+ ooB(t>-WO-ooBBr_l_ vj> ( + 1 V -W) 
^ u-iT+W-sinW IW *ab\T±, U / 

-W+ra^- 1 ^-^)]' ( 6 - 22 ) 
(^-^^Lq)^ = (^L )(^A)^(+l,t;) + (5t-^)(^A)fi a{) (+l^) 

l—^(-l,v-W)} (6.23) 



Ty+2u-27r 



and 

{Ay)-\±L )d z h ah = (±L )(d z A)V ab (+l,v) + (5t-5z)(d z A)n ab (+l,v) 

-[A + 7e - Ci W - ln(2it)] fl^+l, u) + (Si W) $ a6 (+l, v) 

fi a6 (-l, 2u + u) + E(w) $ o6 (-l, 2m + u) - u-^abi-l 

+ ^TW^w i^W K ^(+1, « - WO 

fa^-M-WO], (6. 



M-TT+VK-smPV L W 
2u-2TT+W-sinW 



W+2u-2iv 



where 



A = ln[(47rLo 1 ) 2 Q + ]. (6.25) 

The expressions listed here for h ab and its derivatives neglect additional terms 
of order 5x a ln[(47rLQ 1 ) 2 Q + ] or higher, which are zero in the limit Sx a — > 0, pro- 
vided of course that this limit is taken from directions outside the tangent plane 
to the world sheet at X a (u,v). However, the expressions also contain terms of 
order m[(47rZ/Q x ) 2 (5+] and Sx a d b \ia[(4:7rL^ 1 ) 2 Q + ], which are undefined in the limit 
as 5x a — > 0. The presence of the singular term 4 / uln[(47rLQ 1 ) 2 Q + ] 1 l , ab (+l, v ) in 
h ab is a reflection of the conical singularity that is expected the mark the location 
of any zero-thickness cosmic string, as has been explained in more detail in Paper 
I. Fortunately, none of the singular terms makes a contribution to the first-order 
equation of motion, as will been seen shortly. 

Another important point to note is that equations (16.161) and fl6.2ip - fl6.24p are 
valid only near points X a (u, v) with u G (0, n), which lie on one of the two smooth 
segments of the loop shown in Figure 1. The metric perturbations near points on 
the second segment, for which u G (— 7r, 0), can be calculated by making use of 
the identities I+(u — n,v) = I-{u,v) and — 7r, v) = I + (u,v) for u G (0,7r), 
and the analogous identities satisfied by the pairs (S + , <9_) and (C+, C_). In view 
of equation (I5.37P for h ab , it follows that all the components of h ab except h tz , h xz 
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and h yz are half-range periodic in u near the world sheet (that is, h a b{u — n, v) = 
h a b(u, v)), while the components h tz , h xz and h yz are half-range anti-periodic (that 
is, h a b(u — ii, v) = —h a b(u,v)). This behavior was previously mentioned of the 
tangential projections h a bX^,A (which are all half-range periodic in u) 

in Paper I. Although Paper I deals with a loop that is strictly stationary, the 
statements it makes about the local behavior of h a b are still applicable here, as 
the local effects of loop evaporation first appear at order /i 2 in g ab . 

6.2. The equations of motion 

It is now possible to write down and verify or solve the equations of motion 
(12.131) and (12.141) . Since the only non-zero component of the 2-metric j AB is 
7™ = (j^L )~ 2 (to leading order in /i), the zeroth-order equation of motion 
cfej CD X& ) ^,cD= reads simply q d X^, uv = 0. It is clear by inspection of (14.161) 
that X? y uv = 0, and so (12.131) is automatically satisfied. The second term in 
the first-order equation of motion (12.141) . which is proportional to r y CD 'XjL,cD-> is 
consequently also zero. 

The details of the calculation of the remaining terms in (12.141) can be found 
in Appendix C. Because all the terms in the first-order equation of motion are 
projected onto q d , and so lie in the subspace orthogonal to the tangent plane 
at X a (u, v), the equation contains only two independent components. In what 
follows, it proves to be convenient to decompose each term in the directions of 
the vectors 

M a (s, v) = [1, - sin(f ) x+cos(w) y+s z] a and N a (v) = [0, cos(v) x+sin(w) y] a , 

(6.26) 

which span the orthogonal subspace. 

The first and third terms in (I2.14p . which depend on the perturbation SX a 
and therefore on the correction functions F, G and A, are: 

q c dl CD 5X d , CD = 2fi(±L )-\F uv + sA uv + G u - K )M c + 2 f x(±L o y 1 G uv N c (6.27) 

and 

- 2q^ AC 7 BD Vab X? 0) , A 5X\ B X( 0) , C n = 2^L )- 1 [ks(\u\ - |tt) - (F u + sA u )]N c . 

(6.28) 

The remaining terms depend on h a b and its derivatives, which were evaluated 
earlier in the case u G (0, 7r) only. For positive values of u (or equivalently, for 
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s = +1) it turns out that 



-q c d l Ac l BD h ah Xl 0))A X\ 0))B Xf 0) , CD = 16/i(iL )-i[l n ( W + 2u - 2vr)-ln(2 M )]iV c , 

(6.29) 

-q^j^X^AX^Bh^ = 16 /U (^L )- 1 S(n)[-(2cos2 M )M c + (sin2 M )iV c ] 

+ 16/j(±L )- 1 K(u) [(2 sin 2u)M c + (cos 2u)N c ] 
+M^L o y \_ n J_ smW [(l - cosW)M< - (sin W)N<] 

,a (1_T \-! W+2u-2n-2smW 1 

x[3(l - cosW)M c - (sinW)N c ] (6.30) 

and 



Iq^X^X^Bha^ = -16 yU (^L )- 1 [K(n)sin2n-S( M )cos2n]M c 

g / 1 T \ — 1 W+2u— 2-7T— sinV^ 1— cos W jitc 

^Ita^i M-7r+VK-sinVK V7+2«-2tt 

^(^Lq)- 1 wwjT W o N c . (6.31) 

"^47T u ' u — TT+W — smWW + 2u — 2TT \ ' 

Now, it was shown in Paper I that h a bX^,A X^,b is half-range periodic in u, 
and it is easily seen that taking the inner product with q^) AC 'j BL } X^, C d preserves 
this symmetry. Hence, the functional dependence of the term on the left of (16.291) 
for u G (—77,0) is the half-range periodic extension of the term on the right. 
Similarly, the sum of the terms in (16.301) and (16.311) - which is to say, the term 
on the right-hand side of the first-order equation of motion (12.141) - is just 2a a , 
where a a is the acceleration vector introduced in Paper I. It was demonstrated in 
Paper I that the t-, x- and ^-components of a a are half-range periodic in u, while 
the z-component is half-range anti-periodic, and this is also true of the vectors 
M c (s,v) and N c (v). Hence, the first-order equation of motion requires the sum 
of (16.271) and (16.281) to have the same symmetries, and since the term /€s(|w| — ~7r) 
in (16.281) is half-range periodic this means that F + sA and G must be half-range 
periodic as well. 

These properties can also be predicted from the geometry of the problem. 
The two segments of the evaporating loop should share the same geometry, so the 
correction functions F, G and A are expected to have the same symmetries as the 
corresponding components of the original ACO loop. That is, F and G should be 
half-range periodic in u as the t and x-y components of the ACO loop are, while 
A should be half-range anti-periodic as the ^-component \u\ — \a is. 
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On setting s = +1 in the first two terms (16.271) and (16.281) . the equation of 
motion for u G (0, it) reduces to two partial differential equations: 

F uv + A uv + G u - k = 8[K(u) sin 2u - E(u) cos 2u] + A u _^_ ^ W ^~-t W 

(6.32) 

and 

G uv + k{u - |tt) - (F u + A u ) + 8 [In (W + 2u - 2tt) - ln(2w)] 
= 8[Jf(«) cos 2u + E(«) sin 2u] - 4 ^^^ (sin - 2^^ cos W) 

(6.33) 

which appear as the coefficients of M c and iV c respectively. 
At a schematic level, these equations read: 

F uv + A uv + G u = P^u) and G uv - F u - A u = P 2 (u), (6.34) 

where the functions Pi and P 2 represent all the remaining terms. Taking the v 
derivative of each equation in turn and adding or subtracting the other leads to 
two forced simple harmonic equations: 

G uvv + G u = P^u) and (F u + A u ) vv + F u + A u = -P 2 (u), (6.35) 

whose general solutions are: 

G = ci{p) + a ( u ) cosv + b{u) sin v + J P\{u) du (6.36) 

and 

F + A = c 2 (f ) + b{u) cos-u — a{u) sint> — J P 2 iu) du, (6.37) 

where C\ and c 2 are arbitrary continuous 27r-periodic functions of v, and a and b 
are the restrictions to (0, tt) of arbitrary half-range periodic C 1 functions of u. 

The physical significance of the arbitrary homogeneous contributions to ( 16.361) 
and ( I6.37P can be understood by referring to the intrinsic 2-metric 'Jab- The 
components of 7 AB are expanded to order // in Appendix D.l, where it is shown 
that 7 AB takes on the self-similar form (14.221) at this order if and only if 

F U v + Kv = 0, F uv + F vv + A VV = and G v + F vv = 0. (6.38) 

The first of these constraints reduces the first equation in (16.341) to G u = -Pi (it), 
and so from (I6.36j) a and b are both zero. The trigonometric functions of v in 
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(16.361) and (I6.37P therefore represent free vibrations of the loop with amplitudes 
of order /z, which break the self-similarity of the evaporating loop and, if allowed 
to remain, would contribute to the metric tensor g ab at order /z 2 . 

The second of the constraints (I6.38j) now reduces to F uv = — c? 2 (v), and since 
F is by assumption continuous and half-range periodic in u it follows that c' 2 = 
and therefore that c 2 is a constant (and can be absorbed into J P2(u)du). 
Consequently F uv = 0, and F has the generic form F\{u) + F 2 (v) for some choice 
of continuous and suitably periodic functions F% and F 2 . Furthermore, the third 
constraint now becomes G v = —F 2 , and so from (16.361) c\{y ) = —F 2 . 

The possibility of further constraining F, G and A arises from a consideration 
of transformations of the spacetime coordinates x a = [t, x, y, z] a of the form x a — > 
x a where 

x a = x a + ny a {x b ). (6.39) 

The only transformations of this type that are admissible are those that preserve 
the form of the metric tensor g ab = r] ab + \ih a \, at order //. Since 

9ab = Vab + ^ h ab + ^Vca d bV C + Wcb^ (6.40) 

at linear order in //, admissible transformations satisfy the condition t] c ^ a d^y c = 
and so are generators of the Poincare transformations (that is, combinations of 
boosts, rotations and translations). 

Under such a transformation, the position vector X a of the loop is replaced 

by 

X a = Xfa + 8X a + iiY a (u,v) (6.41) 
where X? s and 5X a are given by (jClO and (10.3j> respectively, and 

Y a (u, v) = ^L [-5F( Uj v), 8G(u, v){cos(v) x + sin(v) y} + 5A(u, v) z} a (6.42) 

for some choice of functions 5F, 5G and SA^ 

It is only possible to impose the admissibility condition r] c ^ a d b )y c = on Y a if 
the derivatives are taken tangent to the worldsheet of the loop, and so Y a must 
satisfy the equations 

X(o), u -Y,u = 0, X (0 ),„ -Y, v = and X {0) , u -Y, v +X {0) , V -Y, u = 0. (6.43) 

6 Note that the x- and y-components of Y a have the forms shown so as to preserve the 
isotropic gauge condition discussed in Section 4, following equation (|4.20l) . 
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In terms of SF, SG and 5 A these equations read, for u G (0, 7r), 

SF u + SA u = SG + SF v = and SF V + SA V + SF U = 0. (6.44) 

The choice SF = —F 2 (v), SG = F^iv) and SA = F 2 {v) satisfies all three equations, 
and specifies a spacetime transformation which reduces F and A to functions of 
u only, and sets c\ = in (I6.36p . 

In a suitable set of spacetime coordinates, therefore, the correction functions 
become 

G = Jp 1 (u)du and F + A = - J P 2 {u) du, (6.45) 

with F a continuous and half-range periodic function of u alone. In fact, F is 
effectively an arbitrary function of u, as it is always possible to define a transfor- 
mation u — ► u' = u + \x\J (u) of the gauge coordinate u, where U is continuous and 
half-range periodic but otherwise arbitrary, under which X a — > X a + fiUX^, u 
and so F —>■ F — U and A — > A + U for it e (0, tt), but G remains unchanged. 



6.3. Solving for F, G and A 

The equation (16.91) specifying u — n as a function of W can be used to eliminate 
u from the terms in (I6.32j) and (16.331) which mix both u and W . The resulting 
expressions for Pi and P 2 read: 

P^u) = K + 8[K(M)sin2M-S(M)cos2M]+4(l-cosI^-W r sinW)/J(I^) (6.46) 

and 

P 2 {u) = -/t(M-|7r)-8[ln{2^ 1 (l-cos^)}-ln(2M)] 

+8[K(u) cos 2u + E(u) sin 2u] - 4(sin - W cos W)/J(W), 

(6.47) 

where 

J(W) = W-\l - cos W) + \W - sin W (6.48) 

and now 

K(u) = Ci{2iy- 1 (l-cosW / )}-Ci(2M) and E(u) = Si{2iy- 1 (l-cos W)}-S 

(6.49) 
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Furthermore, the derivative of (16. 9p reads du = —W 1 J(W)dW, and so 

G = ku + 8 J[K(u)sm2u-Z{u)cos2u)du-4: J W' 1 ^ - cosW - W sin W) dW 
= Ku-A[K(u)cos2u + E(u)sm2u]+4,ln(l-cosW)-4,hi{uW)+C 1 (6.50) 

and 

F + A = \ku(u-ti) + % J[\n{2W~ l {l-cosW)}-\n(2u)]du 

-8 J[K(u) cos2m + £(m) sin 2m] du - A J W _1 (sin W - W cos W) dW 
= \ku{u--k)- %W- l [l + In W - ln(l - cos W)\ (1 - cos - 4W(2 - In W) 
— 8(m Inn — u) — A[K{u) sin 2m — £(m) cos 2m] — 4 / ln(l — cosuj) du + C2. 

J2tt 

(6.51) 

Here, the integral Jln(l — coscj) da; unfortunately cannot be expressed in closed 
form as a combination of standard functions. (The lower bound in the integral 
has been set to uj = 2n because W — > 2n~ as u — > + .) Equations (I6.50P and 
(16.511) can of course be verified directly by differentiation.. However, in tracing 
out the original calculations - that is, in integrating Pi and P 2 - it proves to 
be most efficient to integrate the terms K{u) sin 2m — S(m)cos2m in (16.501) and 
ln(l — cos ly) and S(m) sin 2m + K(u) cos 2m in (I6.5ip by parts with respect to u 
first, then integrate the remainders with respect to W . 

Turning to the initial conditions satisfied by G and F + A, it should be recalled 
that G is continuous and half-range periodic in u and so G(0) = G(tt). Similarly, 
F is half-range periodic in u and so F(0) = F(n). By contrast, A is continuous 
and half-range an^-periodic in m, and therefore A(0) = —A(tt). (Note that the 
combination F + sA, which appears in the leading terms (16 .2 7ft and (16.281) in the 
first-order equation of motion, is half-range periodic for u G [— n, 71"] as required. 
However, it is not necessarily continuous at u = 0.) 

Now, the u — > + lower limit of equation (I6.5ip reads: 

F(0) +A(0) = -167r + 87rln27r + C 2 (6.52) 
while the u — > tt~ upper limit reads: 

F(ir) + A(ir) = 8tt - 8tt 1h(2tt) - 4 Si(2vr) + C 2 . (6.53) 
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(Note in particular that J 2n ln(l — costu) doo = 2nln2 at the upper limit.) 
Subtracting (I6.53j) from (I6.52p therefore gives: 

2,4(0) = A(0) - A(n) = -24tt + 16vr In 2vr + 4 Si(27r) (6.54) 

and so 

A(0) = -127T + 87rln27r + 2Si(27r) « 11.328. (6.55) 
Similarly, addition of ( 16. 52ft and (16. 53ft gives: 

2F(0) = F(0) + F(tt) = -8tt - 4 Si(2vr) + 2C 2 (6.56) 

and so 

F(0) = -4vr - 2 Si(2vr) + C 2 « -15. 403 + C 2 . (6.57) 

The constant of integration C 2 in (16.571) represents a vestige of coordinate 
freedom in the specification of the correction function F that, in view of (16.551) 
is not present in the function A. In fact, different values of C 2 correspond to 
the same physical configuration of the evaporating loop, as the value of C 2 can 
be changed at will by making a suitable Poincare transformation (16.411) . This 
follows from the fact that the choice 5G = 5 A = and 5F any arbitrary constant 
(representing a translation in the coordinate t at order jj) automatically satisfies 
the Poincare constraints (I6.44p . One natural gauge choice for F is to set F(0) = 0, 
in which case C 2 = 47r + 2 Si(27r). 

Consider now equation ( I6.50p for G. Taking the lower limit u — > + (or 
equivalently W — > 2n~) in this equation gives: 

G(0) = C u (6.58) 

while taking the upper limit u — > n~ (equivalently W — » + ) gives: 

G(7r) = 7TK + 4 Ci(27r) - 4 In 2vr - 4 7e + C x . (6.59) 
So it follows that G(0) = G(n) if and only if 

k = [4 In 2tt + 4 7e - 4 Ci(27r)]/7r w 3. 103 7. (6.60) 

That is, iff k — 7°/(4vr) with 7 w 39.0025 the radiative efficiency of the ACO 
loop, as expected. 

Unlike C2, the integration constant C\ in ( I6.50P does have an invariant physical 
significance, and it cannot be arbitrarily changed by making a Poincare or gauge 
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transformation without affecting the form of F. In fact, different values of C\ 
correspond to different values for the length of the evaporating loop as measured 
along curves of constant scale factor e - K M u +^) ; as j s explained in more detail in 
Appendix D.2. If the length L cs f of the loop along such a curve is assumed to be 
independent of // (with L cs f = -^L along the curve u + v = 0) then 

d = -28 - 67 E + 10 In 2tt + 6 Ci(2vr) + Att~ 1 Si(2vr) + Att~ 1 Si(4vr) « -9. 514 3. 

(6.61) 

With C\ assigned this value, the correction function G, which is the half-range 
periodic extension of the function on the right-hand side of (16.501) . is plotted for 
u G [— 7r, tt] in Figure 5a. The half-range periodic extension to [— tt, tt] of the 
function F + A defined by equation (I6.5ip with C 2 = 47T + 2 Si(27r) is plotted 
against u in Figure 5b. As mentioned at the end of Section 6.2, the functions F 
and A are not separately uniquely determined, as different choices of F correspond 
to different choices of the gauge coordinate u. However, the requirement that 
A(0) = — A(tt) = A (with A = —12ix + 8tt In 2tt + 2 Si(27r)) does place a restriction 
on the choice of F. One natural gauge choice is to set A proportional to the 
z-component \u\ — ~7r of the unperturbed position function. Then 

A(u) = -^(\u\-^) (6.62) 

is continuous and half-range anti-periodic on [— tt, tt] with A(0) = A, while from 
(16.51 1) F is the half-range periodic extension to [-7T, tt] of 

F{u) = !««(«- 7T) + f (n-|7r)-8H/" 1 [l + lniy-ln(l-cosiy)](l-cosiy) 
-4W(2 - \nW) - 8(u\nu - u) - A[K(u) sin2u - S(w) cos2n] 

pW 

-4 / ln(l -cosw)^ + 47r + 2Si(27r). (6.63) 

J2-K 

The correction functions A and F in this gauge are plotted against u in Figures 
6a and 6b respectively. 

To summarize, the position function (14.201) with F, G and A suitable exten- 
sions of (16.631) . (I6.50p and (I6.62|) is a solution of the first-order string equation 
of motion ( 12.141) that reduces to the ACO solution in the limit as fi — > 0. Fur- 
thermore, if the integration constant C\ is chosen to take the value ( 16.61!) then 
the solution is the unique self-similar loop trajectory (up to gauge and Poincare 
transformations) with the property that the length of the loop retains its ACO 
value -j^Lq along the curve u + v = (which to leading order in \i corresponds to 
the spacelike slice t — 0). 
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7. Matching the Solution to the Minkowski Vacuum 



By prescription, the evaporating ACO loop radiates away the last of its energy at 
t = ti, the moment when it shrinks to a point at the origin. It is therefore to be 
expected that the metric induced by the loop reduces to the Minkowski vacuum 
at all points in the chronological future of the evaporation point x a = [t^, 0] a . If 
the evaporation point is denoted by pl, the boundary of the chronological future 
I + (pi) is the future light cone Fl of pl, plus pi itself. The weak-field metric 
will match smoothly onto Minkowski spacetime in I + {pl) if an appropriate set of 
junction conditions is satisfied on Fl. I will show explicitly that these junction 
conditions are indeed satisfied later in this section. (Note, however, that the weak- 
field metric remains singular in every neighborhood of pi, and cannot be matched 
onto Minkowski spacetime at this one point.) 

7.1. The surface t — t L = |x| is null to order jj, 

Before proceeding, it is important to first establish that the surface Fl is, to 
order //, just the surface t — tL = |x|, which is of course null with respect to the 
unperturbed metric tensor r] ab . To show that the surface t — tL = |x| is also null 
with respect to the first-order metric tensor r\ ah + hob, it is sufficient to show that 
h a t,(t, x) — > as the field point x a = [t,x] a approaches the surface. 

Suppose then that e is a number in (0, 1), and consider the set of field points 
[£, x] lying to the past of the surface t — ir, = |x| and satisfying the conditions 

\-x\-t + t L = et L and \{t - t L ) 2 + |x| 2 ] 1/2 > s l/2 t L . (7.1) 

The first condition implies that e is a dimensionless measure of the distance from 
the point [t, x] to the surface t — tL = |x| in the background Minkowski manifold. 
In fact, it is easily seen that the Euclidean distance from the point to the surface 
is just -j^Eth- The second condition ensures that the field point avoids a certain 
neighborhood of the evaporation point [t^, 0]. (The scaling factor e 1 ! 2 is chosen 
here because the exclusion of a neighborhood about [t^, 0] with the smaller radius 
et L does not lead to a useful bound on the components of h ab , as will be seen 
shortly.) It is evident that as e — > the points satisfying (17.11) include all points 
in some timelike past neighborhood of every point on the surface t — tL = |x| 
except the evaporation point. 

Eliminating t from the conditions (17.11) yields one non-trivial bound on |x| if 
e < 1, namely 

|x| > l[ £ + (2e~e 2 ) 1 / 2 ]t L , (7.2) 
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which in turn implies that |x| > ttjS £,l ( a s can readily be verified by dividing 
the expression on the right of (17. 2p by e l / 2 ti and showing that the minimum value 
of the resulting function is I/a/2). 

Now, the corresponding image point [t, x] on the future light cone F of the 
origin has coordinates 

t = e K ^|x| and x = e K ^R(-ip)x (7.3) 

where 

e ^ = [( | x |_ f + tj .)/ fi .]-i = e -i (7.4) 
and R is the rotation matrix defined in (15.141) . So in particular 

t = e- 1 |x| > ^e- 1/2 t Ll (7.5) 

while of course |x| = i. 

Note here that i — > oo as e — > 0. Also, according to (15.221) . 

h ab (t, x) = R a c (tfj) h cd (t, x) R b d (tfj) (7.6) 

where the components of the rotation operators R% are all bounded functions of 
ip. Hence, to prove that h a b(t,x) — > as e — > in the conditions (17.11) . it is 
sufficient to show that \h a b{i, x)| is bounded above on Fq by a function of i which 
goes to zero (uniformly with respect to x) as i — > oo. It should also be clear that 
a proof of the last statement would not be sufficient if the factor e 1 ' 2 in the second 
condition in (17. ip were replaced by e or a multiple of e. For then (17.21) would be 
replaced by an inequality of the form |x| > ce for some constant c, and the image 
point [i, x] would be placed under the much weaker constraint i > c, which does 
not guarantee that h a b(i, x) — » as e — > 0. 

Now, given that the components of ^ a b in the equation (I5.36P for h a b(t,5t) all 
satisfy ^^(s, u)| < 1, it follows that 

/■Vl+27T nVo 

\hab(t,x)\ < 4/i / [T+ (Z+ \n) -vj-Uv + ifi / [T - (Z + \n) - v}- 1 dv 

JVo JVi 

= 4 M |/ + |+4 M |/_|, (7.7) 

where [T, X, Y, Z] = AttLq x] as before. So \h a b(t, x)| will be bounded above by 
a function which goes uniformly to zero as t — > oo if the same is true, separately, 
of |/ + | and |/_|. 



38 



Consider then 



rV 

\I_\= [E^v)]- 1 dv where E(s, v) = T - s(Z + ±tt) - v. (7.8) 

JVi 

Since E is positive semi-definite (see Appendix A. 2), there are two possible cases 
that need to be considered: either (i) E(l,v) = for some v G [Vi, Vo], or (ii) 
E(l,v) > for all v G [Vi,V ]. 

Case (i) has been treated in some detail in Appendix A. 2, where it is shown 
that E(l, v ) = for some v only if the point (X, Y) lies on the unit circle (so 
that T = (1 + Z 2 f /2 ) and Z > -\n. Furthermore, if Z > |vr (and so T > 
(1 + jK 2 ) 1 ! 2 m 1.8621, which will henceforth be assumed) then the values of Vq 
and V\ coincide, and according to (1A.9j) 

|/_|=ln[l + 7r/(Z-l7r)]. (7.9) 

In view of the inequality ln(l + x) < x it follows that 

|/_| < n/(Z - |tt) = tt/[(T 2 - If' 2 - |tt] (7.10) 

where the function of T appearing on the right is bounded above by 8tt/T if 
T > ^(4tt + 8V167T 2 + 63) w 2.0869. (The significance of the factor 8?r will 
become apparent shortly.) So, in this case, 

\I-\<8n/T if T> 2.0869. (7.11) 

In case (ii), E(l,v) > for all v G [VI, Vo]. Since E is a linear decreasing 
function of v, there exists a number V* in [Vi, Vq] with the property that E(l,v) > 
\T for all v G (Vi, K) and < E(l,v) < \T for all v G (K, V o )0 An immediate 
consequence is that 

\I-\= [E(l,v)}~ 1 dv+ [E^v)]' 1 dv (7.12) 
JVi Jv, 

where 

v* 

[E(l, v)]- 1 dv < 2(K - Vi)/T < 4tc/T, (7.13) 



7 The value of V* will of course depend on [T, X, Y, Z]. Also, V* may coincide with either Vo 
or V\ , in which case one of the two domains specified for v will be empty. 



v 
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as K - V x < V - V 1 < 2tt. 

The integral over the remaining segment (14, Vo) can be bounded above by 
using equation (I5.34p to replace v with u as the integration variable, so that 

rV r u ( v *) 

/ [E(l,v)]~ 1 dv= / [Hiv^du (7.14) 
Jv* Jo 

where H{v) = T — (u + v ) + Xsint> — Y cosv is a second positive semi-definite 
function, examined in some detail in Appendix A. 3. In particular, it is evident 
from the first line of ( 1A.10I) that 

H{v) > \Z-(\u\ - |tt)| >Z-\tt (7.15) 

if u is in [0, 7r] and Z > ^ir. Also, the condition E(l,v) < \T corresponds to the 
inequality 

Z>\T -\-k-v>\T -^-(l + Ivr 2 ) 1 / 2 (7.16) 

where v here has been replaced by its maximum possible value, ix + (1 + ^n 2 ) 1 ^ 2 
[see Appendix A.l]. In particular, the condition Z > |7T will be satisfied if T > 
16.2906. Hence, if the latter condition is assumed, 

f^* W)]- 1 du < u(V*) [\T - 2tt - (1 + ivr 2 ) 1 / 2 ]- 1 (7.17) 
Jo 

where the maximum possible value of u(V*) is 7T. The function of T appearing on 
the right of (17.171) is therefore bounded above by An/T if T > 8n + 2^4 + 7r 2 w 
32.5811. 

In summary, the inequalities ( 17. 13ft and ( 1A.10I) together imply that 

\I-\<8n/T if T> 32.5811 (7.18) 

in case (ii), and so in view of AT. 11 j) |/-| is bounded above by a function that 
goes to zero uniformly as T — ► oo, as required. Note that the factor of 87T in the 
bounding function arises as a consequence of choosing V* to separate the domains 
E ^ |T. The factor | is to some extent arbitrary here, as any factor less than 1 
will induce similar bounds on |/_|. However, a factor 1 or greater is not viable, 
as it generates no useful bound on Z, and hence no useful bound on the integral 

(EUD. 

The proof that the second contribution |/ + | to \h a b(t,5c)\ is also uniformly 
bounded above by 87r/T for sufficiently large values of T is formally identical to 
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the proof just given, provided that |/_| is everywhere replaced with |J + |, E(l,v) 
with E(—l, v), Z with — Z, Vq with V\ + 2tt and V\ with Vo- Also, the limit u = 
in the integral on the right of (I7.14p should there and below be replaced with 
u = —Ti, and the bound u(V*) < it with u(V*) < 0. 

In conclusion, therefore, \h a b(i, x)| < 64n^/T = 16L /i for sufficiently large 
values of i, and so h ab (£, x) — > as e — > in the conditions ( 17. ip . It follows that 
the future light cone F L of the final evaporation point is, to order /i, just the 
surface t — = |x|. 



7.2. Checking the junction conditions 

The junction conditions that apply when matching two spacetime regions across 
a null hypersurface M were first usefully formulated by Barrabes and Israel [20] 
and have been handily summarized by Poisson in [21]. The first step in the 
Barrabes-Israel method is to parametrize the null hypersurface in the form x a = 
9\9 2 ), where r is an arbitrary coordinate along the generators of the hyper- 
surface, and each generator is labeled by the two remaining coordinates 6 A (for 
A = 1,2). This parametrization can in turn be used to construct an orthonull 
basis {k a , e", e|, iV a } at each point on Af, with 

k a = dy a /dr, e a A = dy a /d6 A (7.19) 

and the remaining (null) basis field iV a uniquely determined by the constraints 

g ab N a k b = l and g ab N a e b A = 0. (7.20) 

The hypersurface then has an associated mass density p, current density j A 
and isotropic pressure p given by 

p=-i<J AB [C AB ], j A = ^a AB [C 3B ] and p = -±[C w ], (7.21) 

where a A B = 9abe A e B and C a p = e£e^Vb-/V a , with a, (3 running from 1 to 3 
and el = k a . Also, [C a/3 ] denotes C Q/3 (A/' + ) - C a(3 (N~) where M~ and M + are 
the (timelike) past- and future-facing sides of Af respectively. (Note that C a p 
typically has different values on the two sides of Af because the limiting values 
of the covariant derivatives V b N a are different, even though the metric tensor g ab 
is assumed here to be continuous across Af.) The two spacetime regions can be 
matched across Af without invoking an extraneous source of stress-energy on Af 
if p, j A and p are all zero. 
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In the case of the evaporating ACO loop, Af is the null hypersurface F L defined 
by t — tL = |x|, and g ab = Vab there. An obvious choice of parametrization is 

y a (r, 6\ 6 2 ) = [t L , 0] a + [r, r sin 6 1 cos 6 2 , r sin 6 l sin 6 2 , r cos 9 1 } a . (7.22) 

The corresponding orthonull basis fields are: 

k a = [1, sin 9 1 cos 6 2 , sin 6 1 sin 6 2 , cos 6 1 ] 11 , (7.23) 

e\ = [0, r cos 1 cos 9 2 , r cos ^ x sin 9 2 , -r sin (7.24) 

e a 2 = [0, -r sin 6 l sin r sin 9 1 cos 0] a (7.25) 

and 

N a = i[l,-sin0 1 cos0 2 ,-sin0 1 sin0 2 ,-cos^ 1 ] a - (7-26) 

In particular, the line element a ab d6 A d6 B intrinsic to the hypersurface is just 
that of a spatial 2-sphere: 

a AB = -r 2 diag (1, sin 2 9 1 ) and a AB = -r~ 2 diag (1, esc 2 6 l ). (7.27) 

To linear order in the derivatives of h a b, the transverse curvature of Fl is 

C aP = e a a e b p d b N a - \e a a e b pN c {d b h ac + d a h bc - d c h ab ), (7.28) 

where the first term e^e b pd b N a is identical on the two sides of Fl, and the remaining 
terms are of course zero on the Minkowski side F£. Because h ab = on F£ the 
derivatives e^d a h bc tangent to F£ are all zero, but the limiting values of the 
transverse derivatives N c d c h ab are not necessarily zero there. So 

[Cap] = -YJ^dJi^-. (7.29) 

Now, the derivatives of h ab can be evaluated formally by using ( 15. 22ft to write: 

e M*,*) = [{itK)KS,*)R b q +RZKS,*) (^R b q )]d c ^+(dx r /dx c )R a p d r h pq (i,5c) 

(7.30) 

where d c i[) and dx r /dx c can in turn be calculated from (15.121) and (15.131) . However, 
as is shown in Appendix A. 3, the derivatives d r h pq (i, x) are divergent at image 
points that lie in the beaming directions of the two kinks, and the loci of such 
points extend to future null infinity on Fq. Although points of this type have 
z = ±g£o, and so (since z = e~ K ^z) the only points on Fl which are limit points 
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of the singularities of d c h a f, lie on the plane z = 0, this feature of the metric makes 
(17.301) less than ideal for calculating [C a p]- 

Instead, it is more convenient when calculating [C a p\ to use the formal defini- 
tion of the derivative 

Jim \\Sy\\- 1 [f(y a + Sy a )-f(y a )-Sy c dJ}=o, (7.3i) 

where / denotes any of the relevant combinations of the components of h a b, and 
1 1 • 1 1 is the Euclidean norm. This approach has the advantage of demonstrating 
that the relevant combinations of the metric derivatives do in fact exist at all 
points on the null boundary (except at pl itself). 

Consider then a fixed point y a = [tL+r, rr] a on F L \{p L } (with f = (x, y, z) and 
|f | = 1), and let x a = y a + 5y a be a field point at a Euclidean distance wt^ from 
(and in the timelike past of) y a , so that 5y a = zut^St, 5r] a where Sr = (Sx, Sy, Sz) 
and \ \[5t, 5r]\\ = 1. If e denotes e~ K ^ (as in the previous section) then 

e = |t L Vf + zu5r\ -Tj}r -wbi= (r • 5r - 8t)w + 0{w 2 ). (7.32) 

It is easily seen that e < \^2zu, and that e > when x a lies in the past of y a . 
The image point x a corresponding to x a has components 

i = £ _1 |rf + zut L 5r \, x = e _1 [(rx + zut L 5x) cost/; + (ry + zut L 5y) sin?/>], 

y = £~^\iry + w ti,5y) cos ip — (rx 4~ w t^bx) sin ip\ and z = e~ x {rz + wt^z). 

(7.33) 

If the non-dimensionalized components [T, X, Y, Z] = 47tLq x, y, z] are substi- 
tuted into the equation (I5.33P for Vk then, after squaring both sides and rearrang- 
ing, the equation reads (for k = or 1) 

\rt+wt L 6t\(V k +hr)-rr-Xk = mt L {5v-X k + ±Kfi(E/tu) [(Vfc + /c7r) 2 -(l + 7r 2 /4)]}, 

(7.34) 

where the factor of k/x arises because j^Lq = KfitL, and 

X fe = (cos{V k + VO, sin(\4 + i>), {k - |)tt). (7.35) 

Equation (17.341) suggests that, for small values of w, the root Vk can be ap- 
proximated by the solution of the equation Vk + kiv — f • which I will denote 
by V k °°. In explicit form, the equations for V^°°,i read 

V™ = x cos(K 00 + V) + y sin^o 00 + V) - p (7.36) 
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and 

V 1 °° + n = x cos(V£° + ip) + y sm^ 00 + ip) + p. (7.37) 

Note however that V^°°,i are typically still functions of w, as ip — —(K[i)~ l \n.E 
increases without bound, and in particular lim ro ^ Vk° does not exist unless x = 
y = 0. 

Although lim ro ^ Vk° does not generally exist, it can be shown that lim ro ^ |^4 — 
V k °°\ = 0, a fact that will prove important in establishing that certain combina- 
tions of the metric components are differentiable on F^. To show that lim ro ^ |^4 — 
V k °°\ = 0, first rearrange ( 17.34p to give 

r(V k + kn - r ■ X fc ) = zut L {5r ■ X k + ^(e/zu) [(V k + kn) 2 - (1 + n 2 /A)]} 

+ (r - \rr + wt L 5r \){V k + kit) 
= (7.38) 

Since | V k + kir\ is bounded above (see Appendix A.l), it is clear that the function 
Q(tu) on the right-hand side of this equation tends to zero as zu — ► 0. The root 
V k °° satisfies the same equation (IT.38j) with 6 (to) = 0. 
Subtracting the equation for V k °° from 07. 38!) gives: 

V k -VF-x[cm(V h +$)-caa(VF+1,)]-y^ 

(7.39) 

In view of the identities cos 2A — cos 2B = 2 sin(Z? — A) sin(A + B) and sin 2A — 
sin 2B = 2 sin(A — B) cos{A + B) this equation becomes 

^-Vr+2[£sin(iy fe +±y^ 

(7.40) 

and so 

\V k - V k °°\ - 2{x 2 + y 2 y/ 2 \ sin[i04 - V?>)}\ < ^Ml. (7.41) 
Since (x 2 +y 2 ) 1 ^ 2 < 1 and \w\ — 2\ sin(|w)| > ^ | w| 3 for all real w, it follows that 

^-^l 3 ^ 1 !©^)! (7.42) 

and therefore lim ro ^ tyk — ^4°°! = as claimed. Furthermore, because Q(zu) goes 
to zero linearly in zu, \V k — V k °°\ is at worst of order zu 1 / 3 for small values of zu. 

The asymptotic forms of the components of h a b can now be developed by first 
substituting the expressions 

X± = ^L Q l e-\t±z)±l (7.43) 
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into dSISHD-dSHOD to give 

h = lL eY^ + O{e% I = ±L e^ + 0(e 2 ), 
S ± = ±±L e cosV l- c z osV » +0(e 2 ) and C± = ±^L e ^f Vl + 0(e 2 ) 

(7.44) 

(assuming for the moment that t ± z ^ 0). 

With t — \rr + zxjt L 5r\ and z = rz + wt L 5z, equation (15.371) generates the 
following expansion for the components of h ab : 

M*>*) = iW) Ml*l,*)itf(v) 

AVz- (l-z)7r (co-ci)f (so-si)z -AV - (1 - z)n 

(cq — Ci)z AV z — {1 — z)it Ci — c 

(s — S\)z AV z — (1 — i)7T Sx — s 

— A^/ — (1 — Z)n Ci — Co Si — So AV"Z— (1 — 5)7T 

(7.45) 



2L fie 
7rr(l — z 2 ) 

+0(w 2 ), 



where 

c fe = cos(T4 + V), s fc = sin(Vfc + V) and AF = - V . (7.46) 
In terms of the components of f the relevant orthonull basis vectors are 
k a = [1, x, y, z] a , el = r(l-z 2 )- 1 / 2 [0, xz, yz, -{x 2 +f)\ and e a 2 = r [0, - ,y. r. 0]" . 

and so 



(7.47) 



k a k b h ab = -4(L /r)fj£ + 0(n7 2 ), k a e\h ab = 0{w 2 ) 
k a e\h ab = ^ (1 J 2)3/2 Lo/ig(Vi + vr - f ■ Xi - V + f ■ X ) + 0(w 2 ) 



(7.48) 
(7.49) 



and 



a AB e a A e b B h, 



ab = -r- 2 [e$e\h ab + (l-f)-h a 2 e b 2 h ab ) 

= -^F) (L /r)fj£[(yx + 7T — f-Xi — Vo + f- X )z - tt(1 - z 2 )] 
+0(w 2 ). (7.50) 

Now, given that lim ro ^ \Vk — V k °°\ = and Ck and S& are continuous functions 
of V k , 

lim(T4+7r-f-Xi) = lim [V^-Hr-ac cos(K 1 °°+^)-y sin(y i °°+^)-f z] = 0, (7.51) 
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and similarly lim ro ^ (K) — f ' X ) = 0. Hence, in view of the expansion (17.321) . 
lim w- l k a k b h ab = -4(L /r)/i(f • 5r - 5t), lim w~ x k a e\h ah = 
lim w~ l k a e b 2 h ab = and lim w' 1 a AB e a A e b B h ab = 4(L /r)/i(f ■ 5r - 5t). 

(7.52) 

Since 5y a = wt L [5t,5r] a it follows from fl7T3TD that k a k b h ab , k a e b A h ab and 
u AB e A e b B h ab are all differentiable at each point on the null boundary F£ with 
z 2 7^ 1, and that 

r A; b ,9 c /i ab = -a AB e a A e b B d c h ab = 4(L /r)^V[l, -*]c and k a e b A d c h^ = 

(7.53) 

there. Given that N c [l, — f] c = 1, it seems therefore that 

[C 33 ] = -^ B [C^] = -2(Lo/r)*2V and a AB [C 3B ] = 0, (7.54) 

and hence that matching the evaporating ACO spacetime to the Minkowski space- 
time across F L induces a boundary layer with a non-zero density p and isotropic 
pressure p. 

However, it should be recalled that L^t^ 1 = Attkli, and so the apparent density 
and pressure of this boundary layer, 

p = —p = —np 2 /r, (7.55) 

are quadratic in Thus p and p are zero at the level of the first-order approxi- 
mation used throughout this paper, and the evaporating ACO spacetime can be 
matched to the Minkowski spacetime without inducing a boundary layer at this 
level of approximation. 

All that remains now is to show that the same conclusion holds at all points 
on F£ with z 2 = 1 (that is, on the z-axis). In this case, x = y = and tends 
to the limiting value 

V k °° = -kn + (k— 1)ttz. (7.56) 

Also, it is convenient to use a slightly different asymptotic form for the spatial 
components of the field point x a , namely 

x = -^L x*e, y = j^L y*e and z = rz + j^L z*e, (7.57) 

incidentally, it is to be expected that neither p nor p depend separately on Lq or i^, as the 
values of the latter will vary as the zero point of t changes, but the stress-energy content of Fl 
should remain invariant under time translations. 
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where x*, y* and z* are all of order e°. 
Then 

X = x* cosip+y* smip, Y = y* cos ip — x* sin if) and Z = 47tLq £ rz+z*, 

(7.58) ' 

and the analogue of (I7.43P becomes 

X± = 47rrL 1 (l±5) £ - 1 +^5(l±5)±i7r+^Lor- 1 (X 2 +F 2 ) £ -|( i LL ) 2 r-^^(X 2 +F 2 ) £ 2 

(7.59) 

to order e 2 , while the analogue of equation (17.381) for can be used to calculate 
the asymptotic expansion 

X±-V k = 4vrrL 1 (l ± z)e~ 1 + z*z(l ± z) ± |tt + A;vr - (fc - |)tt£ 

+ ^L r- 1 (X 2 + F 2 + 2F I + l)e + I(^L ) 2 r - 2 X (2Y + z)e 2 
+ l(±-L ) 2 r- 2 z[(k - |)tt - z*](X 2 + F 2 + 2F z + l)e 2 (7.60) 

(plus terms of order e 3 ). 

The expansions need to be developed to such a high order because if z = 1 
(resp. z = —1) then the terms appearing at order e~ x and e° make no contribution 
to 7_, or C_ (resp. 7 + , 5+ or C + ), and the terms of order e plus the terms 
proportional to (k — ~) at order e 2 dominate the asymptotic behavior.. Otherwise 
(if z = —1, resp. z = 1), it is the terms at order e" 1 and e° that dominate. A 
straightforward but lengthy calculation shows that 

/+ = I- = -\{Lo/r)e, S + = -l(l-z)(L /r)e, 

S_ = \{l + z){L /r)e and C + = C_ = (7.61) 

to order e, and hence from (15.371) that 

+ 0(s 2 ). (7.62) 

ab 

Furthermore, when z 2 = 1 the spherical coordinates (r, 6 |1 ,^ 2 ) parametriz- 
ing the null surface Fl are singular, as <tab is degenerate. A more suitable 
parametrization of Fl in this case would align the poles of the spherical coor- 
dinate system with (say) the x-axis. Then 

k a = [1, 0, 0, z) a , e a x = r[0, 0, -1, 0] and = r[0, 5, 0, 0] a , (7.63) 



h ab (t,x) = (L /r)ne 



—2 z cos ^ zsin^> 

zcos-?/' —2 — cos'0 

5 sinip —2 — sin^ 

— cos if) —sin?/' —2 



47 



and a AB = — r 2 5 B . With this choice of orthonull basis, it is easily checked that 
k a k»h ab = -4(L /r)fie + O(e 2 ), k a e\h ab = 0(e 2 ), k a e b 2 h ab = 0(e 2 ) (7.64) 
and 

a AB e a A e b B h ab = A(L /r)^e + 0(e 2 ),. (7.65) 

just as in the case z 2 ^ 1 considered earlier. So it can be seen that the null bound- 
ary Fl\{pl\ separating the evaporating ACO spacetime from the relic Minkowski 
spacetime contains no stress-energy to first order in /x. 

8. Dynamical Stability of the ACO Loop 

Realistic cosmic string loops would not, of course, have the exact shape of an ACO 
loop, and would not evaporate in a strictly self-similar manner. It is therefore 
important to be able to estimate the impact that deviations from the ACO loop 
trajectory would have on the evaporation process. A full analysis of the stability 
of the ACO loop in the presence of back-reaction is beyond the scope of this paper 
(not least because the dynamical consequences of the self-gravity of a string loop 
in the strong- field limit are not yet understood), and in this section I will limit 
myself to some observations about the dynamical stability of the flat-space ACO 
solution. 

The trajectory of the flat-space ACO loop is fully described by the position 
vector X(r, a) defined in (13.21) . where the mode functions a and b are specified 
by equations (13. lip and (13.121) . For the purposes of the perturbation analysis of 
this section I will refer to the unperturbed mode functions as a and b , so that 

aoK) = - \L) z for - \L < a + < \L (8.1) 

and 

b (a_) = ^[cos(27ra_/L) x + sin(27ra_/L) y], (8.2) 

where a± = r ± a. (Strictly speaking, of course, a is the even periodic extension 
of the function on the right of (18.11) .) 

A generic perturbation of the ACO solution has mode functions of the form 

a(cr + ) = a (cx + ) + <5a(cr + ) and b(a_) = b (cr_) + 5b(cr_), (8.3) 

where 5& and Sh are continuous functions of their arguments, and are both small, 
in the sense that |5a| L and |<5b| L everywhere. Without loss of generality it 
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can be assumed that the perturbed solution has the same parametric period, L, as 
the unperturbed ACO loop, and so 5a and 5b are periodic functions with period L. 
The only other constraints on 5a and 5h are the gauge conditions |a'| 2 = |b'| 2 = 1, 
which to linear order in the perturbations require that a' • 5a' = and b' ■ 5h' = 
at all points where these derivatives exist. 

Note in particular that the derivatives 5a' and 5b' need not be continuous, nor 
is it necessarily true that either \5a'\ or \5h'\ is small. A simple way to visualize 
the effect of a perturbation is to plot the unit vectors a' and b' on the surface 
of a unit sphere, an approach first popularized by Kibble and Turok [22] . Figure 
6 contains three plots of this type. The first represents the unperturbed ACO 
loop, and shows the b' curve following the equator of the sphere while a' = ±z 
is concentrated at just two points, the north and south poles. In the second plot, 
the perturbations 5a' and 5b' are assumed to be both continuous and small, with 
the result that the a' and b' curves depart only minimally from the unperturbed 
curves. The a' curve, in particular, is understood to jump discontinuously from 
the small circle around the south pole to the small circle around the north pole 
when cr + « then back again when cr + « tt, so the two kinks are preserved under 
perturbations of this form. 

The third plot in Figure 6 shows the effect of a perturbation in which 5a' is 
continuous but no longer small. In fact, |5a'| is assumed to be of order unity 
for a narrow range of values of the parameter <r + centered on a + ~ and a + »s 
it, creating two smooth segments on the a' curve linking the north and south 
poles. This has the effect of removing the discontinuity in the a' curve, and thus 
smoothing the two kinks, but an unavoidable consequence is that the a' and b' 
curves now cross at two points, indicating that a pair of cusps will momentarily 
form during each oscillation of the perturbed loop. 

So although the unperturbed ACO loop does not support cusps, it is evident 
that there exist solutions of the flat-space equation of motion arbitrarily close 
to the ACO loop in a-b trajectory space that do support cusps. The presence 
of cusps complicates the back-reaction problem considerably. Because the bulk 
velocity of a string loop momentarily reaches the speed of light c at a cusp, a non- 
negligible fraction of the loop's total energy can be concentrated there. A simple 
kinematic argument ([23], Section 6.2) suggests that the total energy M r inside 
a sphere of radius r centered on a cusp scales as r 1 / 2 , and so the gravitational 
potential M r /r diverges as r -1 / 2 . This in turn is an indication that the effects of 
gravitational back-reaction at a cusp cannot be modeled accurately at the level 
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of the weak-field approximation]^] 

Any cusps forming on a perturbed solution that is close to the ACO loop will 
necessarily be very narrow, as the energy in such a cusp will be proportional to 
the range of the parameter cr + covered by the segment linking the two poles on 
the Kibble- Turok sphere. What effect back-reaction would have on the cusps is 
still an open question. It is possible that back-reaction would act to narrow the 
cusps even further, thus driving the perturbed solution closer to the unperturbed 
ACO loop. The ACO loop would then be dynamically stable. 

Alternatively, the cusps could broaden in the presence of back-reaction and 
ultimately expand to the size of ordinary macrocusps, which typically contain a 
fraction of order /x of the total energy of the string loop ([23], Section 6.2). In this 
case, it seems likely that the perturbed loop will evaporate very quickly, as cusps 
are known to beam away copious amounts of gravitational radiation [19J. Indeed, 
the presumption alluded to in Section 3 that the ACO solution has the lowest 
radiative efficiency of any cosmic string loop encourages the view that dynamical 
instabilities in the trajectory of the ACO loop, if they exist, are likely to be of 
only secondary importance. A perturbed trajectory that diverges rapidly from 
the ACO loop will have an increased radiative efficiency and will evaporate more 
quickly, whereas perturbations that remain close to the ACO loop will survive 
longer. 

A third possibility is that cusps on perturbed solutions close to the ACO loop 
are unstable to gravitational collapse, with consequences that have not yet been 
explored and can at present only be imagined. 

Another potential source of dynamical instability is loop self-intersection. 
String loops that intersect themselves are believed to split into two daughter 
loops with high probability ([23], Section 2.7). The ACO loop itself has no self- 
intersections, but there do exist trajectories arbitrarily close to the ACO loop that 
intersect themselves, as will be seen shortly. 

A self-intersection occurs on a general flat-space loop when 

a(a+ + A) - a(<r+) = b(<r_) - b(a_ - A) (8.4) 

for some values of cr+. <r_ and A, with A G (0, L/2) ([23], Section 3.8). If a and 
b are decomposed in accordance with (18.31) . the x- and ^-components of equation 

9 By contrast, because the bulk Lorentz factor of a string is bounded at a kink, the energy 
M r inside a sphere of radius r centered on a kink typically scales as r, and the gravitational 
potential M r /r has no problematic divergences. The breakdown of the weak-field approximation 
at a kink is due not to a local spike in the energy content of the string, but to the fact that the 
tangent vector X? Q -., U is discontinuous (see Section 9 below). 
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QgID read: 



cos(2ttct_/L) - cos[2tt(ct_ - A) / L] = 2nL-\5 A a - 5 A b) ■ x (8.5) 

and 

sin(27nr_/L) - sin[27r(<r_ - A)/L] = 2nL- 1 (5 A a-5 A b) ■ y (8.6) 
respectively, where 

5 A a = 5a(a + + A) - 5a(a + ) and 5 A b = 5b(o~_) — <5b(o~_ — A). (8.7) 
Equations (I8.5P and (18.61) can be solved to give: 

cos(2vra_/L) = nL~\5 A a - 5 A b) ■ [x+^g^y] (8.8) 

and 

sin(27nx_/L) = ttL" 1 ^ - 5 A b) ■ [y - *] ■ (8.9) 

In particular, 

1 - cos(2ttA/L) = 27i 2 L- 2 [\(5 A a-5 A b) ■ x| 2 + \{5 A a - 5 A b) ■ y| 2 ], (8.10) 

and so A must be of order \5 A a — 5 A b\ or smaller. For small perturbations of the 
ACO loop, therefore, self-intersections can only occur with A C L. 
Moreover, for small values of A equations (I8.5P and (I8.6P reduce to 

A^ 1 (5 A a— 5 A b)-x. w — sin(27rcr_/L) and A~ 1 (5 A a—5 A b)-ymcos(2iT(7-/L) 

(8.H) 

respectively, and can be satisfied only if either 5a' or 5b' is discontinuous or of 
order unity at the self-intersection point (cr + , <r_) ~ (cr + + A, cr_ — A). Hence, at 
least one of the a' and b' curves will deviate significantly from the corresponding 
unperturbed ACO curve on the Kibble- Turok sphere. 

The picture that emerges, therefore, is that a flat-space solution close to the 
ACO loop can intersect itself only by pinching off a small daughter loop with 
parametric period A C L. A solution that does just this is easily constructed 
by choosing (at time t = say) a small segment of the unperturbed ACO loop 
with length of order A and deforming it so that it crosses itself. But the same 
perturbation procedure could be applied to any flat-space loop solution, so self- 
intersections of this type are in a sense rather trivial. 

Furthermore, the larger of the two daughter loops created by the self-intersection 
will be little changed from the original parent loop. It will have a marginally 
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smaller parametric length L — A, and in place of the self-intersection there will 
be two narrow kinks traveling in opposite directions around the loop. Analytic 
studies of kinks on long strings strongly suggest that gravitational back-reaction 
will suppress small-scale structure on a string with a characteristic length r on a 
timescale At ~ rj (87r 2 /z) ([23], Section 6.7), and this estimate has been confirmed 
in numerical simulations of planar loops carrying up to 32 kinks by Quashnock 
and Spergel [II] . In the case of a kink formed by the pinching off of a daughter 
loop with a parametric period A the characteristic length r will be of order A, and 
the decay time At for the kink will be much smaller than the radiative lifetime 
t L = L/ (47T/t/i) of the ACO loop. So it is to be expected that the two narrow kinks 
created by the self-intersection will quickly dissipate, leaving behind a daughter 
loop that very closely resembles the ACO loop. 

To summarize, there is as yet no reason to believe that the ACO solution is 
dynamically unstable, although a definitive statement on the matter will depend 
on the results of full back-reaction stability analyses, both at the weak-field level 
(without cusps) and in the strong- field limit (with cusps). Nonetheless, it seems 
safe to conclude that the effect of possible self-intersections on the stability of the 
ACO loop is negligible. Furthermore, even if it were to turn out that the ACO 
loop is unstable to certain mechanisms (such as cusp formation), this would not 
in itself invalidate the importance of the self-similar evaporating ACO solution 
described above, as the longest-lived loops in any ensemble of loops with the 
same initial energy would be those that most closely resemble the ACO loop - on 
the presumption, of course, that the ACO solution is unique in having the lowest 
radiative efficiency of any string loop. 

I should add that at present there is an ongoing debate about the nature of 
the scaling peaks in the production spectrum of cosmic string loops, fueled by 
the conflicting results of recent numerical [2U ES] and analytic [2S1 EZ] studies of 
the primordial string network. However, this debate has no direct bearing on the 
stability or evolution of the ACO loop. 

9. Conclusions 

In this paper I have constructed a solution of the linearized Einstein and Nambu- 
Goto equations describing a cosmic string loop that evaporates self-similarly under 
the action of its own self-gravity, eventually radiating away all its energy and 
momentum to leave a remnant flat spacetime, and whose spacelike cross-sections 
are identical (to leading order in the string's mass per unit length /i) to those of 
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the Allen-Casper-Ottewill loop. 

From a purely physical viewpoint, the ACO loop is one of the most important 
flat-space cosmic string solutions, as it is long-lived and approximates well the 
lowest-efficiency daughter loops observed in low-resolution string network simu- 
lations [IE1 I19j 10 l . The fact that its evolution and ultimate evaporation are an- 
alytically tractable is an added bonus, and makes the ACO loop an object of 
potentially great interest to mathematical relativity. 

However, much work needs to be done before it could be claimed with any 
confidence that the evaporating ACO loop can be described by a self-consistent 
solution of the Einstein equations. A first step in this direction would be to 
construct a solution to the strong-field back-reaction problem that reduces to the 
bare ACO loop in the limit \x — > - or, at the very least, to demonstrate that 
such a solution exists. A second desideratum would be to find a solution of the 
Abelian Higgs field equations that can act as a suitable source for the evaporating 
ACO metric. Neither task is likely to be easy. 

One question that arises naturally in this context is whether the self-similarity 
of the linearized solution would be preserved in the corresponding fully non-linear 
solution (if indeed such a solution exists). This question remains an open one, but 
it can be given a more rigorous formulation as follows. If the original Minkowski 
coordinates x a = [t,x] a describing the linearized solution are written in terms of 
the similarity coordinates x], where if) and x are defined by equations (15.121) 
and (I5.13P respectively, then it can be shown that the vector field 

k a = dx a /dif) = [0, —y, x, 0] a - Kfj,[t - t L , x, y, z] a (9.1) 

satisfies the equation 

(Vc(a + h c (a)d b) k c + \k c d c h ab = -n^r] ab (9.2) 

10 More recent high-resolution simulations indicate that string loops, when they first form in 
an evolved string network, have an approximately fractal structure, with a fractal dimension 
approaching 2 on large scales [281 EH EQ]. The ACO loop therefore approximates realistic 
incipient string loops less closely than was first believed. Nonetheless, it remains true that the 
ACO loop has the lowest known radiative efficiency of any string loop. Furthermore, as was 
explained in Section 8, there are good theoretical reasons for believing that radiative back- 
reaction acts to preferentially eliminate small-scale structure from a loop on timescales much 
shorter than the lifetime of the loop (although the principal decay channel for structure on the 
very smallest scales would probably be particle production [IE])- Hence, it seems likely that 
any fractal loop similar in shape to an ACO loop on large scales would radiate away much of 
its small-scale structure and quickly come to resemble the ACO loop on all scales. It has also 
been suggested [26l [27] that loop fragmentation will quickly act to suppress fractal structure, 
by triggering a cascade of daughter loops with progressively fewer cusps and kinks. 
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to linear order in fi. That is, k a satisfies the linearized version of the equation 
V( a fcb) = —Kfig a b for a conformal Killing vector field. 

The corresponding non-linear solution will therefore preserve self-similarity if 
it admits a vector field k a satisfying the exact conformal Killing equation V( a A;&) = 
—K\xg a b-, in which case it is always possible to find coordinates [V>,x] in terms of 
which the metric tensor has the conformal form 

^(^,x) = e- 2 ^V(x). (9.3) 

(The conformal Killing field in this case is k a = 8^.) So the question of whether 
self-similarity would be preserved by a fully non-linear solution can be recast as a 
question about the existence of solutions to the strong-field back-reaction problem 
with the conformal form (19.31) . 

Another of the more problematic aspects of the ACO loop is the behavior of the 
gravitational field near the two kink points. The discontinuity in the null tangent 
vector X^, u at the kinks generates two apparent anomalies in the weak- field 
analysis described in Section 6. The first is that the term X^,cd appearing in 
the first-order equation of motion (I6.28P is undefined at the kinks, and so it cannot 
meaningfully be said that the equation of motion is satisfied at these points. The 
second anomaly is that the metric perturbations h a b fail to be differentiable on a 
manifold of points extending from the trajectories of the kinks out to future null 
infinity, as is explained in more detail in Appendix A. 3. 

The breakdown in the differentiability of h a b poses no intractable difficulty 
for the linearized analysis of Section 6, because the ill-defined terms make no 
overall contribution to the linearized Einstein tensor. However, the fully non- 
linear Einstein tensor contains products of terms of the form d c g a b, and in the 
generic case it is expected that a distributional interpretation of the Einstein 
equation is possible only if the metric derivatives d c g a b are locally square integrable 
[2]. So the breakdown of differentiability in the beaming directions of the kinks 
may bode ill for the prospects of a strong-field solution. 

A further problem is that the non-differentiability of h a b is superimposed on the 
conical singularity marking the string worldsheet at the kink points themselves. 
The singularity in h a t, at the kinks was examined briefly in Paper I, and is far 
more pathological than the simple logarithmic singularity that appears at all 
other points on the string worldsheet. An appreciation for the nature of the 
singularity can be gained by repeating the analysis performed earlier, in Section 
6.2, to calculate the near-field metric perturbations at a kink point. 
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If the field point x a = -^L [T, X,Y, Z] a is chosen to be close to the lower of 
the kink points (at u — 0) near the spacelike surface t = 0, so that 

[T,X,Y,Z] a = [v,cosv } sinv, -ln} a + AttLq^x" 1 (9.4) 

for some v E (— it, ir), then it can be shown that x± = v + 47tLq x (5t ± Sz), and 
that 

X+ - V x - 2tt = ^(4ttL 1 ) 2 g_, x- - V x = 2ir - BkL^Sz 

and x±-V = 4xL?(6t ±8z-\ ^f^J (9.5) 
to first or leading order in £x a , where 

Q_ = (5x) 2 + (5y) 2 + (6t + 5z)(5t + Sz + 25a; sin w - 25ycosv) 

= y ab (-l,v)5x a 5x b . (9.6) 

All four of the functions in (19. 5p appear as arguments of the logarithmic and 
sine and cosine integral functions contributing to the terms I±, S± and C± in 
(15.371) . For points on the future lightcone of the kink, where r] ab Sx a Sx b = 0, the 
expansions (19.51) simplify considerably, as x± ~ then becomes 47tLq x (5t ± 5z), 
while Q_ = 2(5t + 5z)(5t + <5xsint> — 8ycosv). (Furthermore, the function 5t + 
Sx sin v — 5y cos v appearing in Q_ is zero if and only if the field point x a + 5x a 
lies in the beaming direction of the kink: see Appendix A. 3.) 

However, for the purposes of calculating the near-field derivatives d c h ab it is 
necessary to retain the term r] ab 8x a 8x b in full, with the result that many of the 
derivatives are undefined in the limit as 5x a — > 0. The near- field expansion of 
h ab appearing in the footnote in Section 4.2 of Paper I is as complicated as it is 
because the "off-shell" contributions have been included so as to give a hint of 
the complexity of the derivatives. 

It should be clear from the foregoing remarks that the breakdown of the field 
equations at the kinks is not strictly speaking a strong- field effect, but stems 
rather from the zero-thickness assumption underpinning the Nambu-Goto action, 
which imposes infinite worldsheet curvature at the kinks. It seems unlikely that 
the singularity at the kinks would be resolved in a fully non-linear solution of 
the Einstein equations, and I suspect that the kinks can be smoothed only by the 
inclusion of a finite-thickness field-theoretic source. However, the fact - mentioned 
in the footnote in Section 8 - that the Newtonian gravitational potential M r /r 
remains bounded near the kinks (and in fact is likely to be of order fi there) offers 
some hope that the kinks on the ACO loop could be incorporated into a fully 
non-linear solution. 
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There are other topics of interest in this area that would repay further investi- 
gation. For example, it would be useful to be able to offer a definitive proof that 
the ACO loop does indeed have the lowest radiative efficiency of all flat-space 
cosmic string loop solutions, or to demonstrate that the evaporating ACO loop is 
at least neutrally stable to (cusp-free) small perturbations. 

I would like to thank the referee for highlighting the need for some helpful 
clarifications. 



A. Some Properties of the Roots V\ and Vq 
A.l. Bounds on V\ and Vq 

If the vector (X, Y, Z) corresponding to the spatial components x = -^L Q (X, Y, Z) 
of the image point is represented in standard spherical polar coordinates: 

(X,Y,Z) = (r sin 6 cos0, r sin 8 sin 0, r cos#), (A.l) 

then since T = r the equation (15.331) for the roots Vk specialized to V\ reads 

Vi + 7r = r — [(r sin 6 cos0 — cos Vi) 2 + (r sin 8 sin — sin Vi) 2 + (r cos# — \^) 2 } 1 ^ 2 ■ 

(A.2) 

Differentiating the expression inside the square root with respect to V\ indi- 
cates that V\ + 7r takes on an extremal value whenever sin(l / 1 — 0) = - un- 
less r sin 9 = 0, in which case equation (1A.2I) solves explicitly to give V\ + 7T = 
r — [1 + (r =p \' n ) 2 ] 1 ^ 2 (with the sign choice =F corresponding to cos 6* = ±1 if 
r / 0). The case rsin6* = is of course equivalent to X = Y = 0. and is 
discussed in more detail at the end of Section 5.2 It is easily verified that, in this 
case, V\ + 7r is a strictly monotonically increasing function of r, and so lies in the 
range [-(1 + ivr 2 ) 1 / 2 , -|vr) if cosfl = +1 and in [-(1 + \n 2 ) 1/2 , |tt) if cos# = -1. 

If sin(Vi — </)) 7^ then the right-hand side of (IA.2I) is extremized if either 
(cos Vi, sin Vi) = (cos 0, sin 0) or (cos Vx, sin V]) = (— cos0, — sin0). It follows 
from fTA~2|) that 



r-[(rsin^ + l) 2 + (rcos^-|vr) 2 ] 1/2 < Vi + tt < r-[(r sin#-l) 2 + (r cos^-iyr) 2 ] 1 / 2 . 

(A.3) 

Again, differentiation with respect to r of the lower and upper bounding functions 
in (IA.30 is sufficient to show that both are monotonically increasing functions of r 
(and furthermore the left-hand function is strictly monotonic unless tan# = — 2/tt, 
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while the right-hand function is strictly monotonic unless tan# = 2/tt). From this 
fact it is readily seen that 

- (1 + K) 1 / 2 < V x + 7T < (1 + ivr 2 ) 1 / 2 (A.4) 

(the lower bound being achieved when r = or tan# = —2/tt, and the upper 
bound when tan 6* = 2/n and r > (1 + j^ 2 ) 1 ^ 2 )- Since the root Vo lies in the 
interval [Vi, V\ + 2tt], it is evident that Vo is also bounded: 

- (1 + \tt 2 ) 1/2 - 7T < V < (1 + |vr 2 ) 1/2 + 7T « 5.00369. (A.5) 

The fact that Vb and VI are bounded should come as no surprise, as V and 
Vi + tt are the values of the scaled time coordinate T at the points of intersection 
of the backwards light cone of the image point [|x|,x] with the trajectories of the 
lower (u = 0) and upper (u = ±7r) kink points respectively. 

A. 2. The zeroes of T — s(Z + \n) - v 

According to (I5.34p . the function v(u) describing the curve of intersection T of 
the world sheet with the backwards light cone of [|x|, x] = -^L [T, X, Y, Z] has a 
possible stationary point whenever E = T — s(Z + ^n) — v = 0. (Recall here that 
s = sgn(w), or equivalently s is the sign of z ■ X, u on the segment of the string loop 
at which the intersection point occurs.) It is important to be able to locate those 
points [|x|,x] on the surface F at which E goes to zero for some value of v in 
[Vi, Vi + 2ir], because E appears as the denominator in the integrals contributing 
to /i a6 (|x|, x) in (15.361) and so the components of h ab are potentially divergent at 
such a point. 

As mentioned in section 5.2, E is positive semi-definite by virtue of equation 
( 15.311) . which can be rewritten in the form 

T — s(Z + |-7r) — v = u - s(Z + |tt) + [(X- cos v) 2 

+ (Y-smv) 2 + {Z-(\u\-±n)} 2 ] 1 / 2 
> u - s(Z + ±tt) + \Z - (\u\ - ±tt)\, (A.6) 

Here, equality occurs if and only if X = cosv and Y = sinv, while the function 
on the third line is zero if and only if u < s(Z + ^n). 

Suppose therefore that the point (X, Y) lies on the unit circle, so that (X, Y) = 
(cosvp, sin vp) for some angle Vp. Then T — s{Z + jir) — vp will be zero whenever 

Z = A — \ix and T = v-p + sX, (A. 7) 
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where A is any number satisfying A > w if w > 0, or A < — u if u < 0. (Note 
that if the null constraint T = (X 2 + Y 2 + Z 2 ) 1 ^ 2 is also imposed then vp = 
[1 + (A — ^tt) 2 ] 1 / 2 — s\, and for each value of s there is a one-parameter family of 
spacetime points on Fq for which E = at some point on T.) 

If < A < 7r then the spacetime point -^Lq[T, X, Y, Z] coincides with the 
point X a (u, v) = X a (s\, vp) on the world sheet. It follows from the results of the 
previous paragraph that E will be zero for all values of u in the range [0, A] if 
s = 1, and all values of u in the range [— tt, —A] if s = — 1. This in turn means 
that the curve T includes a horizontal segment with v = vp over the ranges of u 
indicated. 

More important for the purposes of calculating /i a b(|x|,x) are image points 
that lie off the worldsheet. These have \Z\ > ^n, and so either A > n if s = 1 
or A < if s — — 1. In geometric terms, such points either lie on the upwards 
[Z > |7r) extension of the helical segment of the ACO loop with u > 0, or on the 
downwards [Z < —\n) extension of the helical segment with u < 0. Furthermore, 
it is evident from the preceding analysis that in this case E will be zero on the 
curve T for all values of u G [0, 7r] if s — 1, and all values of u G [— tt, 0] if s — — 1, 
and so T includes a straight-line segment with v = vp over one of these ranges. 
An immediate consequence is that V = Vi = vp if s = 1, and V = V\ + 27r = t> F if 
s — — 1. (This straight-line segment corresponds to the "trivial solution" = v 
mentioned in Section 6.1.) 

It is evident now that if the image point [|x|,x] lies on one of the two helical 
extensions of the loop, but not on the loop itself, then the corresponding integral in 
the equation (15.361) for /z ab (|x|, x) is undefined. However, the problematic integrals 
can be reformulated by making use of the equation (15.341) for dv/du to replace v 
with u as the variable of integration. Furthermore, since (X, Y) = (cos vp, sin vp) 
the denominator T — (u + v ) + X sin v — Y cos v becomes simply T — u — v p, where 
vp = x± = T ± (Z + |7r) if s = =Fl. Hence, 

Vi+2tt 

[T + [Z + jti) — f] _ll I , a fe(— 1, v) dv is replaced by 




[-(Z+Itt)-u) l ^ ab {-l,x+)du 



and 

/ [T — (Z + |7r) — w] _1 \I / a b(-|-l, v ) dv is replaced by 

JVi 



V 
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The new integrals are easily evaluated, giving 
1+ = -]n[(±n-Z)/(-fr-Z)], S + = J + sin( X+ ), C + = J + cos( x+ ) (A.8) 
(if s = — 1) and 

/ = - In[(Z + |tt)/(Z - |tt)], SL = I_ sin( X _), C_ = /_ cos( X _) (A.9) 
(if s = 1) in place of the corresponding terms in fl5.38p - fl5.40p . 

A. 3. The zeroes of T — [u + v) + A sin v — Y cos v 

A second consequence of f)5.34p is that the graph of v against u on the intersection 
curve T typically has a point of vertical inflection whenever H = T — (u + v ) + 
X sinv — Y cosv = 0. As was indicated in Section 5.2, the factor H is positive 
semi-definite, because 

T — [u + v) + Xsin-u — Y cost" = [(AT sin i; — K cost") 2 + (Xcosw + Y sinf — l) 2 

+{Z-(\u\ - ±7i)} 2 } 1/2 + X sinv - Y cosv 
> | A sin f — K cosf | + Asinf — K cosf , (A. 10) 

with equality occurring if and only if Acosf + Y sinv = 1 and Z = \u\ — |7T. 
In particular, it follows that spacetime points with H = somewhere on T have 
\Z\ < |7r. It was seen earlier that spacetime points outside the world sheet with 
E = somewhere on T have \Z\ > |7r, so the two classes of points are disjoint 
outside the world sheet. 

Suppose now that H = at the point (u,v) = (uf,vf) on T. The equations 
H = and Acosf + Y sinv = 1 (with Asinf — Y cosv < 0) can then be solved 
exactly to give 

T = u F + v F + (R 2 - 1) 1/2 , X = cosv F - (R 2 -l) 1/2 sinv F 
and Y = sinv F + (R 2 - l) 1/2 cosv F , (A.ll) 

while of course Z = \u F \ — |7T. The value of R > 1 is fixed by the null constraint 
T = (A 2 + Y 2 + Z 2 ) 1/2 , which reads explicitly 

u F + v F = [R 2 + (\u F \ - ivr) 2 ] 1 / 2 - (R 2 - l) 1 ' 2 . (A.12) 
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So there exists a two-parameter family of spacetime points with the property that 
H = somewhere on T. 

If [T, X, Y, Z] is held constant and u varied from up to up + 5u, a perturbation 
expansion of the equation (15.281) for v indicates that 

v(u F + Su) -up = -(65u) 1/3 (A.13) 

to leading order in Su. Hence, v remains a monotonic decreasing function of u 
even in the singular case where H = somewhere on T. 

The importance of the zeroes of H lies in the fact that the spacetime derivatives 
B c h a b of the metric perturbations h a b(i,5t) are potentially divergent at any image 
point [i, x] = ^L [T, X, Y, Z] with the property that H = when t> = Vo or 
Vi. This can be seen by first differentiating the equation (I5.36P for h ab and then 
integrating by parts to give 

/•Vi+2tt 

d c hab = -4:fi(d c x+) (x+ ~ v)~ l V, ah (-l,v)dv 

fV 

-Mdcx-) / (x- - vy'tiabi+i, V) dv 

JVi 

-4/i ( X+ -Vi- 27r)- 1 ^ afe (-l, V l + 2tt) ^(^ - X+ ) 

+4/i( X _ - v 1 )- 1 ^ ab (+i, V 1 )d c (V 1 - x-) 
-4//(x_ - ^ )- 1 ^ a6 (+l, K )9c(Vo - x_) 
+4/x( X+ - V r o)~ 1 *ot(-1, - X+), (A.14) 

where again x± — T ± (Z + f) and fi a b(s,f) = ^^^(s, f ) as defined in Section 
6.1. 

The two integrals on the right-hand side of ( 1A.14I) converge because all the 



components of f2 Q & are identical to components of ^/ a b, and it was demonstrated 
in the previous sub-section that integrals of this type will converge unless the field 
point lies on the worldsheet itself. Also, the derivatives B c x± = 47tLq [1, ±z] c are 
everywhere well defined. 

However, the terms involving d c V\ and d c Vo in (1A.14I) are potentially problem- 
atic. Implicit differentiation of the equation (I5.28j) defining v gives: 

d c v = [T — (u + v ) + Xsinf — Y cos?;] -1 

x4ttL x [T - (u + v), -(X - cosv), — (Y - sinu), -{Z - (\u\ - \tt)}] c , 

(A.15) 
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where the denominator is of course just H . The derivative B c V k therefore typically 
diverges if H = at v = V k (and in fact the components of the 4- vector multiplying 
H^ 1 i n (IA.15I) will all vanish only if the point -^Lq[T, X, Y, Z] lies on the world 
sheet) 

In fact, if k = or 1 and the image point x a = -^L [T, X,Y, Z] a has the 
property that 

XcosVfc+YsinVfc = 1, T = kir+V k —X sin V k +Y cos V k and Z = (\k\-\)-n 

(A.16) 

(with X sin V k — Y cos V k < 0) then the perturbative solution of the equation 
(E33|) for V k at a nearby point x a + 6x a = ±L [T + 5T,X + 6X, Y + 5Y, Z + 5Z} a 
reads: 

V k (x a + 6x a ) - V k (x a ) = Q 1/3 (5T + 5XsmV k -5YcosV k ) 1/3 

+ \{X sin V k -Y cos V k )- l Q 2/3 (5T + 5X sin V k - 5Y cos K fc ) 2/3 , 

(A.17) 

plus terms linear in Sx a . 

Substitution of the expansion flA. 17|) into the equation 05. 36!) for h ab then gives: 

h tz (x a + 5x a ) - h tz (x a ) = -8zfi(5l /3 + \5 2 k /3 ) + 0(5x), (A.18) 

h xz {x a + 5x a ) - h xz (x a ) = -Sz^sinVk + \6 2 k /3 (smV k + 2D k cosV k )] + 0(5x) 

(A.19) 

and 

h yz {x a + Sx a ) - h yz {x a ) = 8zfx[5l /3 cos V k + ^J /3 (cos V k - 2D k sin V k )} + 0{Sx), 

(A.20) 

where 

5 k = 6D k 3 (5T + 5X sin V k - 5Y cos V k ) , (A.21) 

n It might be supposed from (|A.14[) that d c h a i, is also undefined whenever x± — V = 0, with 
V = Vi, Vq or V—\. That this is not so (at least for image points outside the world sheet) can 
be seen by replacing v with u as the integration variable in the equation (|5.36[) for h a b to give: 

H-^abi-l, v) du - 4/i / H- l ^ ab (+l, v) dv 

-IT JO 

It is clear that d c h a b involves only terms of the form fl a bd c v and H~ 2 d c H 1 where d c H is a linear 
combination of d c T, d c X, d c Y and d c v, and according to (|A.15|) d c v is undefined only when 
H = 0. 
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D k = -(X sin V k - Y cos V k ) > (A.22) 

and z = sgn(Z) = sgn(A; — ~), while all other components of h a b(x a + 5x a ) — h a b{x a ) 
are of linear order or higher in the components of fccf^l 

In summary, if the field point x a satisfies the constraints (1A.16j) for either 



k = or k = 1 then the derivatives of the metric components h tz , h xz and h yz 
with respect to t, x and y are all undefined there. Field points of this type have 
z = ig-^o an d so have the same altitude as one of the kink points on the string. 
The failure of h a b to be different iable at such points therefore has no effect on the 
analysis of the first-order string equation of motion in Section 6, as the equation 
of motion is there developed and solved only for points on the string away from 
the kinks. However, the singularities in the derivatives of h a b do potentially affect 
the matching of the weak-field solution to the Minkowski vacuum across the null 
surface F^. Fortunately, the relevant combinations of the components of h a b are 
all differentiable on Fl, as is shown in detail in Section 7.2. 

It might seem puzzling that the weak-field metric induced by the ACO loop is 
singular on a manifold of points that extends away from the worldsheet. However, 
this pathological behavior has a simple physical explanation. It is well known ([23J, 
Section 6.3) that a generic kink on a zero-thickness string emits a narrow beam 
of gravitational radiation in the direction of its motion. Although the metric 
derivative is singular on the locus of the beam, the singularity is nonetheless 
relatively benign, as the geodesic equations for a free particle remain integrable 
across the beam. 

In the case of the ACO loop, the kink points correspond to u = kn for k = 
or 1 and so trace out the two families of spacetime points 

X ( a 0) = -LL [k7r + v,cosv,smv,(\k\ - \)n\ a . (A.23) 

The instantaneous 4-velocities of the kink points are therefore 

d v Xl Q) = £Lo[l,- sin v, cost;, 0]<\ (A.24) 



12 Incidentally, the factor D^ 3 appearing in the equation for S k is guaranteed to be well defined 
if the field point lies away from the string worldsheet, because if X smVk — Y cos Vk — then 
the previous constraints (|A.16j) imply that 

X = cosV k , r = sinVfc, T = kir + V k and Z = \kir\ - ±tt 

and so the field point is identical to the point (u, v) = (kn, Vk) on the worldsheet. 
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and the beams, which trace out null curves with tangent vector d v X? \, can be 
represented parametrically in the form 

[T, X, Y, Z] = [tor + v, cosv, sinv, (\k\ — ~)7r] + A[l, — sinv, cosv, 0], (A. 25) 

with A an affine parameter. Hence Z = (\k\ — |)7r and, after eliminating A, 

T - (kn + v) + X sinv -Y cosv = 0. (A.26) 

The last equation is just the constraint H = with u = kn. So a spacetime point 
lying off the worldsheet satisfies the constraints (IA.16I) if and only if it lies in the 
beaming direction of a kink. 



B. Calculating h a i and its Derivatives 

With x± = v + (it ± u) + 4ttL 1 (St ± Sz) and 

x+ - V = 2u + SnL^Sz, x+-Vi = W + 2u + AnL \5t + 5z- SV) 

x _ -V = Iu- 1 (AtiL 1 ) 2 Q + and X- ~Vi = W + InL^tft- 5z - SV), 

(B.l) 



where 5 V is of linear order in Sx a and is defined in (I6.13p . while Q + is of quadratic 
order and is defined in (I6.15I) . the expressions (I5.38l) - (l5.40p become, for u > and 
5x a small, 



L 


H = ln(W + 2u- 2tt) - 




7_ = 


= ln[(47rL - X ) 2 Q 


- s + - 




— cos(2w + f ) 


. C +. 




sin(2u + f ) 



_ ] ^ St + Sz-SV Sz 
W + 2u-2ir ~ ~u 

5t-8z- SV 



W 



(B.2) 
(B.3) 



+47TLQ 1 



cos(2u + u) 

sin(2w + v) cos(2w + v) 
cos(2w + t>) — sin(2w + v 



-47rL "?/ 



i„.-i 



suit) 
cost" 



Si (W + 2u - 2vr) - Si (2u) 
Ci (W + 2w - 2tt) - Ci (2m) 

" Si (W + 2u - 2?r) - Si (2u] 
Ci (W + 2u- 2tt) - Ci {2u 

W) 



5z + AttLq 1 (W + 2u-2ti 



sm[v 
cos(f — W) 



{St + Sz) 



(St + 5z- SV) 
(B.4) 
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and 



S- 



{ln[(4 7 rL 1 ) 2 Q + ]+7 E -ln(2 M )} 
+4ttL 1 {ln[(47rL l ) 2 Q + ] + 7e - ln(2u)} 

{St - Sz) 



sint> 




— COSf 


sin v 




~ SiW ~ 


COSf 




sin v 


cosv 




CiW 



cosv 
- sinv 



(St - Sz) 



-AttL, 1 



sin v cos v 
cos v — sin v 



SiW 
CiW 



sin(v — W) 
cos(v — W) 



(St-Sz-SV). 



(B.5) 



Note here that Cix = lnx + 7 E + 0(x 2 ) for small values of x, where 7 E ~ 0.5772 
is Euler's constant, while Six = x + 0(x 3 ). 

If the terms which tend to zero as Sx a — > in (1B.2j) - (lB.5l) are neglected, 
and the resulting leading-order expansions substituted into (I5.37p . the metric 
perturbations h ab take on the form shown in equation (16.16p . 

Also, given that 

d(SV)/d(Sx a ) = (u-n + W -sinW)- 1 

x [u — it + W, {cos(v — W) — cos v } x + {sin(t> — W) — sin v} y + (it 



the terms in the small-distance expansions of the derivatives of I±, S± and C± 
which do not vanish in the limit as Sx a — > can be read off directly from equations 
(E2])-(JE5D: 



u) z\ a 
(B.6) 



8J+ = 47rLn 1 



{W+2u-2ir)(u-ir+W-s\riW) 

dj- = d a A + AttL q 1 



J D a -4vrLo 1 M- 1 [0,z] a 



W(u-n+W-sinW) 



E„ 



(B.7) 
(B.8) 



a S+ 



sin(2w + v) cos(2w + v) 
cos(2w + v) — sin(2-u + v) 



-AttL^u- 1 



sm v 
cosv 



[0, z] a + 4vrL 1 



Si (W + 2u - 2tt) - Si (2u) 
Ci (W + 2u- 2tt) - Ci (2m) 

sin(t> — W) 



i 



(ty+2u-27r)(u-7T+ty-sin W) 



[l,z]a 

cos(v — W) 



D a 
B.9) 
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and 



O a S- 
d a C- 



d a A 



smv 

COS V 



+47rL W ( u _ n+w _ sinW ^ 
where A = \n[(4nL 1 ) 2 Q + ], 



+ 4 7 rL 1 {A + 7 E -ln(2u)} 
{St - Sz) - 4ttL 1 
E n . 



COSf 

— sin v 



[-,--]« 



COSf 

— sinv 



sin v cos v 
cos v — sin v 



SiW 
CiW 



sin(t> — W) 
cos(f — W) 



[l,-z]« 
(B.10) 



D a = [- sin W 7 , — {cos(v— W)— cos-u} x-{sin(f-W / )-sin u} y+ (2m— 27T+W— sin W 7 ) z] a 

(B.ll) 

and 

E a = [sin W, {cos(v — W) — cos i;} x + {sin(v -W) — sin y + (W - sin W) z] a . 

(B.12) 

Again, the expressions fl6.2ip ~ f16.24p for the spacetime derivatives of h a b at 
points on the null surface F and close to the worldsheet of the loop follow directly 
by substituting the leading-order expansions (IB. 71) - (IB. 101) into (15.371) . 

C. Deriving the Equations of Motion 

If the position function X a of the evaporating loop is assumed to have the form 
(gZZQD then 

I ( fl 0) (n 1 .) = lL [« + ^r(« 1 «)r (C.l) 

with 

r(u, v) = cos(f ) x + sin(t>) y + (\u\ — -7r) z, (C.2) 
and the first-order perturbation <5X a is as given in (14.231) . namely 

SX a (u, v) = -±-L Kfi[Uu+v) 2 , {u+v) r(u, v)} a +±L of i[-F, G{cos{v) x+sin(t;) y}+A z] a 



(C.3) 



The first derivatives of X^ are therefore: 



X? 0) , u =±L [l,sz] a and X? 0) , v = ±L Q [1, - sm(v) x + cos(^) y] a , (C.4) 
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with s = sgn(w), while the only non-zero component of X? q \,ab ia^l: 

Xfa, vv = ±L [0,-co S (v)±-sm(v)y] a . (C.5) 
Similarly, the first derivatives of 5X a are: 

5X a , u = -±L k/j,[u+v, r(u, v)+s(u+v) z} a +j^L fi[-F u , G u {cos(v) x+sin(t>) y}+A u 

(C.6) 

5X a , v = -±L K{i[u + v, r(u, v)] a + ±Lofi[-F v , G v {cos(v) x + sin(v) y}+A v z] a 
+ ±L ^{G - k{u + v)}[0, - sin(v) x + cos(t;) y] a (C.7) 

and the only second derivative of 6X a contributing to the first-order equation of 
motion is: 

5X a lUV = -^L K/i[l,sz] a + ^L 0/ u(G u - k)[0, - sin(V) x + cos(V) y] a 

+±L Q fi[-F uv , G uv {cos(v) x + Bin(u) y}+A uv z] a . (C.8) 

The 4-dimensional tangent space at the point x a = X a (u, v) is spanned by the 
vectors X^, u and X^, v , plus 

M a (s, v) = [l,- sm(v) x+cos(w) y+s z} a and N a (v) = [0, cos(v) x+sin(w) yf 

(C.9) 

In calculating the various contributions to the first-order equation of motion 
(12.141) . it is useful to note that all the terms are transvected with q d , where q a b = 
— ^ a b(s, v) is the projection operator orthogonal to the worldsheet at X a (u, v), so 
that 

q c d X d 0) , u = q c d Xf 0) , v = 0, q c d M d = M c and q c d N d = N c . (CIO) 
Other helpful identities involving q d include: 

q c d [0,-sm(v)x+cos(v)y] d = M c , q c d [l,0] d = -M c and q c d [0, s z] d = M c . 

(C.ll) 

13 The claim that X? Q y uu = is true away from the kink points at u = and u = w, but it can 
be seen from the expression for , u that the z-component of X? \ , uu contains a distributional 
term at the kinks. The possibility of pathological behavior at the kinks will be ignored here, 
as the ramifications of the singularities at the kinks have already been discussed in Paper I. It 
seems unlikely that the dynamical role of the kinks can be clarified without developing a solution 
of the Abelian Higgs field equations as a finite-thickness source for the ACO loop. 
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Thus, recalling that the only non-zero component of ^ AB (at zeroth order in 
H) is 7™ = (±L )~ 2 , the first term in f[2TT4j) is: 

q c dl CD 5X d )C D = 2(±L y 2 q c d 6X d }UV 

= 2^L )-\F UV + sA uv + G u -k)M c + 2^L Q )- l G uv N c . 

(C.12) 

The second term in (I2.14p is identically zero by virtue of the zeroth-order equation 
of motion, while the third term is: 

- 2q d >y AC >y BD ri ab X( 0) , A 5X b , B X d 0) , CD = -2( 1 ^LQ)~ A q c d Xf ld) , vv (r] ab Xf Q) , u 5X b , u ) 

= 2 f i(±L r 1 iKs(\u\ - §tt) - (F u + sA u )}N c . 

(C.13) 

The remaining terms in (I2.14p are more complicated, as they involve h ab or 
its derivatives. However, evaluation of these terms is simplified by making use 
of the following table, which represents the inner products of the terms in the 
first column with those in the first row in the form {k\, fc 2 , ^3, ^4), which in turn 
is shorthand for kiX^ a , u +&2^(o)a,i> +k^M a + k±N a . Also, c 2 and s 2 stand in 
for cos 2u and sin 2u respectively, while cw and sw are shorthand for cos W and 
sin W. 





Y b 




^ab(s,v) 










(0,4,-2,0) 





®ab(s,v) 





(-2,0,1,0) 


^ ab (~s,2u + v) 


(-2,0,2,0)c 2 + (0,0,0,-2)s 2 


(-2,-2,3,0)c 2 + (0,0,0,-l) S2 


Q ab (s,v) 





(0,0,0,-1) 


Q ab (-s,2u + v) 


(2,0,-2,0)s 2 + (0,0,0,-2)c 2 


(2,2,-3,0)s 2 + (0,0,0,-l)c 2 


Tlab(s) 





(2,0,-1,0) 


^ab(s) 


(2,4,-2,0) 


(2,2,-3,0) 







(2,0,-l,0)(l-c w ) + (0,0,0,l)s w 


V ab {- S ,v-W) 


(2(l-c w ),4,-2(2-c w ),2s w ) 


(2,2,-3,0)(l- CH /) + (0,0,0,l) Sw 



Furthermore, the vectors X? Q s, u and X? Q ^, V are both null and orthogonal to M a 

and N a , with X {0 ),u -X(p), v = (^o) 2 - 

It should of course be remembered that the expansions (16.161) and (I6.2ip - 
(16.241) were developed for the case u > only, and so s = +1 in the expressions 
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in this table. Because the term involving ln[(47rLQ 1 ) 2 Q+] in the equation for h ab 
is proportional to ^/ ab (+l,v), the term containing /i^X^iX^B in f)2.14p is 
free of logarithmic singularities. Similarly, the singular terms proportional to 
^(+1,^) in fl6.2ip -fl6.24p do not contribute to the term Xt,, B (h bd , a -\h atnd ) 
on the right-hand side of f)2.14p . However, there remain singular terms propor- 
tional to Q ab (+1, v) in the derivatives of h ab whose contributions to the first-order 
equation of motion need to be examined more carefully. 
Now, the fourth term in the equation of motion is: 

~~ Qdl AC l BD habX^,AX^, B X^ }C d = —(-^L ) A q c d XfQy vv (h ab X^, u X^, u ) 

(C.14) 

where the only non-zero term in h ab X^, u X b Q y u is proportional to the coefficient 
oflW-l): 

h ab Xf 0)m X\ 0)m = 16M^L ) 2 [ln (W + 2u- 2tt) - ln(2u)], (C.15) 

and so 

-q c d l AC l BD h ab X^, A X b (0) , B Xf 0) , CD = W^Lo)- 1 ^ (W + 2u - 2vr)-lnf»]iV c . 

(C.16) 

Furthermore, the right-hand expression in (I2.14p can be broken into two parts: 

— g cd >y B Xfo,A X(ty,B h bd , a = —(j^Lq) 2 q cd (X( Q j lU X^, v +X^, V X(ty, n )h bd , a 

(C.17) 

and 

\(f 7 X^,aX^,b h a b,d= ii^Lo) 2 q c X^, u X^, v h ab , d . (C.18) 
Here, q cd X^, u X^, v h bd , a involves terms of the form 

- q cd X b {0)JV d t h bd = 4/i(Si W )M C + A/i(8t - 8z) (d t A)N c + 4/i[A + 7e - Ci W - \n{2u)]N c 

+4 / u£(w) [-(3 cos 2u)M c + (sin 2u)N c ] + AfJC(u) [(3 sin 2u)M c + (cos 2u)N c ] 
+^ u -Jw W sinw M( 1 ~ ^sW)M c - (sinW)N c ] 

- ^iwZinWW + L2, [3(1 " cosH/)M c - (sinH/)iV c ] (C.19) 

and 

- q cd X\ 0))V d z h bd = -4/i(Si W )M C + 4//(<ft - 8z) (d z A)N c - 4/i[A + 7e - Ci W - \n(2u)}N c 

+4 / u£(m) [-(3 cos 2u)M c + (sin 2u)N c ] + A^K(u) [(3 sin 2u)M c + (cos 2u)N c ] 
+^ u ™ + w-lw M( 1 ~ cos W ) MC - (sm^)iV c ] 

+^ 2 7-::w-Zw w + l-2, [3(1 " cosH/)M c - {m*W)N*\. (C.20) 
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So 

— (l^-^o) 2 q cd X(ty, u X(ty, v hbd, a — 1 q cd X h ^, v (d t h hd + d z hbd) 

= A^lL )- l {(5t-5z)[(d t k) + (d z A)]N c 
+2E(u)[-{3cos2u)M c + (sin2u)N c ] 
+2K(u) [(3 sin 2u)M c + (cos 2u)N c ] 
+ u -, + w- smW {( 1 - cosW0M c - (sin 1^] 
+ 2U u-:tw~Jtw W w + L-2, [3(1 - cosW0M c - (sin W)N°}}, 

(C.21) 



where 



; <9 t A) + (d z A) = Q-'dt^^x^x") + Q^d^a^x^x") 

= 2Q~ 1 (y at + y az )5x a = 0, (C.22) 



with here denoting ty a b(+l,v). 

"(0)^^(0) 



Similarly, g cd X f a w„ ATL, U hb d , a is a sum of three terms proportional to 



-q cd X b {0) , u d t h bd = 4/i£(u)[-(2cos2m)M c + (2sin2u)iV c ] + 4fiK(u)[(2sm2u)M c + (2cos2w)iV c ] 

-^ u-,:w W ^ww + l-2, [(2 " cosW)M< - (smW)N% (C.23) 

- q cd X b (0) , u d x h bd = -^ T-^Z-lTw w + l-2, [(2 - cos^)M c - (smW)N<} 

(C.24) 

and 

- g cd X b m , u d y h bd = -8/ C^^ +2 U, [(2 - cosWW* - (smW)N% 

(C.25) 

So 

- (^L o ry d X ( a 0) ^ X (o)>« h bd,a = -{^L )- l q cd X\^ u [d t h bd - (smv)d x h bd + (cosv)dyh bd } 

= 4 / u(^L )~ 1 S(n) [-(2 cos 2u)M c + (2 sin 2u)N c ] 

+M±L )- 1 K(u)[(2sin2u)M c + (2cos2u)N c }. (C.26) 

Adding together (1C.21[) and (1C.26j) gives equation (16.301) . 

The various terms contributing to the second part of the last term in (12.141) 

are: 

X^, u X b 0) , v d t h ab = 8/i( 4^0) [K(u) sin 2u-nu) cos 2u] -8 M (^L ) U _^ W M ifefgfc , 

(C.27) 
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va yb a l _ o,,( 1 r \ cos(u— WQ-cosj l-cos W (n no\ 

^ i0 ),u^( ),vO x tl ab - -S/H^oJ u _„ +w _ BiaW W+2 n-27r' 

and 



X ( a 0) , u X ( 6 0))U ^/i a6 = 8/i(iL )[if(u) sin2w-£( M ) cos 2u\ +8/i( ^L ) 2 :_^ +VK _ sin w w+2m _ 2 ^ ■ 

(C.30) 

Now, if Vd is any contravariant vector and s = +1, then 
q cd V d = -[V t - (sinv)V x + (cosv)V y + V z ] M c - [(cosv)V x + (sinv)V y ] N c . (C.31) 
Hence, the second part of the last term in (I2.14p becomes 
(±L )- 2 q cd X^ u X b (0) , v h ab , d = -IQ^L^Kiu) sin2w - E(u) cos2 M ] M c 

o. I 1 t \— 1 W +2u-2-K-jivaW l-cos i/c 

"^Utt-^O; u -ir+W—BmW W+2u-2tt 

-8u(f Lq)' 1 wwl c ? ° s ^ (C.32) 

D. The Intrinsic 2-Metric and the Length of the Evaporat- 
ing ACO Loop 

D.l. The intrinsic 2-metric 7 AB 

To linear order in /1, the intrinsic 2-metric 7 AB = g a bX a ^X h ,b consists of three 
basic terms: 

lAB — X(o),A -X(p),B +2X( ),(A -5X, B ) +habX^ , A X( ) , B , (D-l) 

with and 5X a as given in (1C.ll) and (1C.3I) . In particular, X^, u ■X^, u = 
X(o), u -X(o), u = and X^, u -X^, v = (^L Q ) 2 , while 

2X {0) , U -5X m = 2fi(±L ) 2 [K(u - ±tt) - (F u + A u )}, (D.2) 

2X( ),(« = ^L ) 2 [k{u - |tt) - 2k(m + u) - F u - + A,)] (D.3) 



2u-27r+H/-sin W l-cos W 



and 

9.Y,^ .AIT = -9.nl 

> 4-7T 

for u > 0. 



2X (0) „ -5X,„ = -2 / u(^Lo) 2 (G + F„) (D.4) 
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One of the components of the term proportional to h ab , namely 

h ab X^, u X b {0) , u = l^{±L Q ) 2 [\n{W + 2u-2n) -ln(2w)], (D.5) 

has already been evaluated in Appendix C. The other two components are: 

h ab Xfo, u X{ 0) , v = -8fi{±L ) 2 {K{u) cos2w+£(«) sin2«-[ln (W + 2u - 2tt)-1h(2u)]} 

(D.6) 

and 

h ab X{ 0) , v X b {0) , v = -8v(±L ) 2 { lE -CiW + K(u)cos2u + X(u)sm2u + lnW 
-[\n(W + 2u-2ir) -\n(2u)]}. (D.7) 

So 

luu = (^Lo) 2 e- 2K ^[-2 / i{F u +A tt -/ t ( M -i7r)}+16/i{ln(^+2 M -27r)-ln(2 M )}], 

(D.8) 

luv = (±-L ) 2 e- 2 ^ u+ ^[l-^{F u + F v + A v -K( U -^)} 

-8//{- \n(W + 2u- 2tt) + ln(2u) + K{u) cos 2u + £(u) sin 2m}] 

(D.9) 

and 

lvv = (^L )V 2 ^)[-2 / i(G + F ?; ) 

_8/i{ 7E - Ci W + In W - \n(W + 2u - 2vr) + ln(2w) + K(u) cos 2u + £(«) sin 2m}] 

(D.10) 

to linear order in u, provided that Kfi\u + v \ <C 1. Here, the factor of q- 2k K u + v ) 
inserted into ^ uv absorbs the term — 2k(u + v) in (ID. 31) . while multiplying j uu and 
7 W by e - 2K v(u+v) com mits an error of order fj 2 . 

Since W is a function of u alone by virtue of (I6.9p . the requirement that 
7^8 take on the self-similar form (j-Lo) 2 e~ 2K ^ w+v '^y AB (u) to linear order in u, as 
explained in Section 4, places three constraints on the correction functions F, G 
and A, namely: 

F U v + A uv = 0, F uv + F vv + A VV =0 and G v + F vv = 0. (D.ll) 

It should be noted that although equations ( ID. 81) - (ID. 101) were derived on the 
assumption that Kfi\u + v\ <C 1, they do in fact hold for all values of u + v. 
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This can be seen by first taking the corrected position vector (14.201) and us- 
ing this to calculate X, A -X,b, with the multiplicative scale factor e ~ 2K f J -( u + v ) re- 
tained but all other terms truncated at order [i. This recovers the contribution of 
X {0) , A -X i0) , B +2X i0) , {A -5X, B) to (|El-([DlOD as shown. 

The contribution of the remaining term h a bX^, A X^,b to ^ AB - where X? ^, A 
is given by (15. ip - can be evaluated by referring back to the discussion of rotating 
self-similarity in Section 5.1. As was demonstrated there, the value of h ab at an 
arbitrary point 

x a = X? 0) (u,v) = ±L [t L (l - e- K ^ u+ ^),e-^ u+v ^(u,v)] a (D.12) 
on the worldsheet is related to its value h ab at the image point 

x a = ±L e K ^- K ^ u+v) [\r{u,v)\,R{-t{j)r{u,v)} a (D.13) 
on the future light cone Fq of the origin through the formula 

h ab = R a c ^)h cd R b d (*P), (D.14) 

where, from (I5.13p . 

tfj = u + v - (/tyu) _1 ln(l + K/i\r(u,v)\). (D.15) 

Here, the image point 

x a = ±L (l + K f i\r(u,v)\)- 1 i\r(u,v)\,R(-^)r(u,v)] a (D.16) 

is just the point X? Ju, v) on the worldsheet with u = u and v = v — ip = 
(k/x) -1 ln(l + Kfi\r(u, v)\) — u, and it was shown in Section 5.2 that k/jl\u +v\ <C 1 
for reasonable values of \i. So 

X^ A (u, v) = X? 0) , A (u, v + ij) = e-^ u+ ^R b %iP) X b (0) , A {u, v) (D.17) 

to leading order in /z, and therefore 

habX(Q} , A ,B I (u,v) =e 2 U+V ^h ab X( }, A X(Qj,B\(u,v)- (D.18) 

That is, to linear order in \i the term h ab X^ , A XL , B at a general point on the 
worldsheet differs from its value at the image point x a only by a multiplicative 
factor of e ~ 2K ^ u+v \ and equations (jn^ - pXTOj) are recovered in full. 
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D.2. The length of the loop as measured along curves of constant scale 
factor 



Equations flD.8j) - flD.10l) can be used to calculate the length of the evaporating 



loop as measured along any curve of constant scale factor e K ^( u + V ) . The length 
L cs{ of the loop along the curve u + v = + v* is given by 

PIT 

L cs{ = 2 (2-f uv --f uu --f vv ) 1/2 \ v=Ut+v ^ u du, (D.19) 



o 



where, to linear order in /i, 

(27™-7™-7,J 1/2 = V2(^L )e-^ u+ ^[l + ^(G + A u ) + 2 f ,{ lE -CiW 

+ \n(2uW) - \n(W + 2u- 2tt) - K{u) cos 2u - sin 2m}]. 

(D.20) 

Hence, in view of the expression (I6.50p for G: 

L csi = 2V2(±L )e- K > x( - u * +v * ) Al + \h(ku + Ci + A u ) 

Jo 

+2m{ 7e + In W - Ci W - 2K(u) cos 2u - 2E(u) sin 2u}} du, 

(D.21) 

where 

J [In W - Ci W - 2K(u) cos 2u - 2E(w) sin 2u] du 

= W~ 1 (2 + \nW)(l - cosW) - ±W(l + ]nW) - \W~\l - cos2W) - § sin W 
-{u - tt) Ci W - Si W + Si(2W) - [K{u) sin 2u - S(u) cos 2u] (D.22) 

and so 

i'JT 

/ [In W - Ci VP - 2K (u) cos 2m - 2S (u) sin 2m] = vr In 2vr + vr - vr Ci (2tt) - Si (4tt) . 
Jo 

(D.23) 

Furthermore, it follows from (16.601) and (16.551) that 

/ Kudu = i/iK7r 2 = 27rln27r + 2tt7 e - 27rCi(27r)] (D.24) 
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and 

/ A u du = A(n) - A(0) = 24vr- 167rln27r-4Si(27r). (D.25) 
Jo 

Substituting flD~24j) and (1TX251) into (jEMJ) therefore gives: 

L csi = (^ 5 L )e~^ (M * +,; * ) [l + |/i{C 1 +28+67 E -101n27r-6 Ci^vr)^^ 1 Si^Tr)^" 1 Si(47r)}]. 

(D.26) 

Note here that L cs f — > -^L Q in the flat-space limit — > 0. This is to be 
expected, as the curves of constant scale factor are curves of constant t in this 
limit, and the speed |<9 t X| is equal to l/v2 at all points on the unperturbed 
ACO loop (except at the kink points at u = and it, where it is undefined). The 
measured length L cs f = -^L is therefore just the Lorentz-contracted fundamental 
length L . 

If /i is non-zero then L cs f will deviate from its ACO value -^L on the curve 
u + v = by a term of order fi unless C\ takes on the particular value 

C x = -28 - 67 E + 10 In 2tt + 6 Ci(27r) + 4tT 1 Si(2vr) + Ait' 1 Si(47r) w -9. 514 3. 

(D.27) 

There is of course no necessary reason to choose this value for C±, but it bears 
repeating that different values of G\ will lead to different "dressed" values of L cs f 
along any given curve of constant u + v. 



D.3. Invariant length of the evaporating loop 

Another measure of the length of the loop that is potentially of interest here is 
the invariant length Li, which is a gauge- invariant quantity defined at each point 
(«*, f*) on the worldsheet by 

L l (u*,v*) = 2( f ^dudv) , (D.28) 

where T>(u*,v*) is the closure of the subset of the worldsheet that is causally 
disconnected from the reference point (u*,v*) [31]. Calculating L\ to order \i 
at a general point on the trajectory (14.201) of the evaporating ACO loop is a 
lengthy process whose details are of marginal importance, so I will offer only an 
abbreviated summary here. 

In the flat-space limit fi = the domain V{u*,v*) is bounded by the null 
curves u = u* and v — v* which, in view of the identification of (u, v) with (u + 
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2irk, v — 2irk) for all integers k, intersect at the points (it*, v* + 2tt) = (it* + 2tt, v *) 
and (u*, f* — 27r) = (u* — 2n, i>*). The boundary of U(u*, v*) is therefore completed 
by the curves 

f = u + 27r — + f * and v = u — 2n — + (D.29) 

and so 

/ / (±L ) 2 dvdu + / / (±L ) 2 dvdu) 

Ju t -2-K Jv* Ju, J U-2lT-U t +V t J 

= L (D.30) 

In order to calculate L\ when fi is non-zero, it is useful to first transform to 
a set of gauge coordinates (V, v') with the property that the curves u' = and 
v' = are (to order ji) null curves through the reference point (u,v) = (w*,t>*). 
A suitable choice is 

u' = u — /i(t> — v*)G{u) and v = v — \iT{u\ u*), (D.31) 

where 

g{u) = G(M)+4{7 E -CiW r + ln^-ln(^ + 2M-2vr)+ln(2M) 
+K(u) cos 2u + £(w) sin 2u) 
= Ku + A( lE -CiW + \nW) + C 1 (D.32) 

and 

r-u 

T{u\u^) = / [{F n + A n -K(M-i7r)}-8{ln(W / + 2M-27r)-ln(2n)}]dn. (D.33) 

In the primed gauge the 2-metric components ( lD.8l) -( JTXT0l) become j u , u , = j v , v , = 
and 

j 1 / 2 ee lu , v , = (±L ) 2 e- 2K ^[l - f,{F u - k(u - ±vr)} + M („ - v*)d u G{u) 
-8/i{- \n(W + 2u- 2tt) + ln(2u) + K(u) cos 2u + E(it) sin 2m}] 

(D.34) 

to linear order in //0 



14 Strictly speaking, the functions of u on the right-hand side of (|D.32|) . (|D.33p and (|D.34|1 
should be the half-range periodic extensions of the expressions shown. However, all three equa- 
tions are accurate if u and both lie in [0,7r]. 
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Recalling again that the point (u, v) is identified with (u ± 2tt,v =p 2n), the 
curves u' = and v' = v* intersect at the points 

[u, v) = (tt* — Am, v* + 2n — Av) = (it* + 2n — Au, — Av) (D.35) 

and 

(u, v) = (u* + An, v* - 2tt + Aw) = (w # - 2vr + Au, v* + Aw), (D.36) 

where 

Am = -2ttu^K) (D.37) 

and 

Av = — uJFfw* + 27r; u*) 

= -2/x / [{F u + A u -K(M-i7r)}-8{ln(iy + 2n-27r)-ln(2u)}]du 

= -[167r-8Si(27r)]u^ -38.92u. (D.38) 

The deviation of v by the amount ±38.92/i from its flat-space values v* and t>* ±27r 
at the four intersection points accounts for the rotational phase shift of the same 
magnitude, mentioned at the end of Section 3, that is induced by gravitational 
back reaction on the ACO loop over the course of a single oscillation period. 

The boundary of V{u*, v*) is completed by the straight line segments (in {u, v) 
coordinates) joining {u* — Au,v* + 2i\ — Av) to (u* — 2n + Am, w t + Av), and 
(w* + 2tt — Au, t>* — Av) to (u* + An, v* — 2n + Av). In the primed gauge, these 
segments become 

v > = v* + (u-u*+2n)[l-n- 1 Av+n- l Au+/ig(u)}+Av-Au-iiJ r (u;u*) (D.39) 
and 

v ' = v* + (u'-u*-2n)ll-n- 1 Av+n- 1 Au+fig(u')]-Av+Au-[iJ : '(u-,u*) (D.40) 

to first order in u. 

The equation (ID.28P therefore becomes: 

1 /2 

pVi(u') /"U*+27T — Au pv* 

Li(w*,0 = 2[/ / ^ 1/2 dv'du'+ / -f 1/2 dv'du' 

(D.41) 



76 



where V\ and v 2 are the functions on the right-hand sides of (ID.39j) and (1D.40P 
respectively, and the form of 7 1 / 2 is given in flD.341) . with u and v everywhere 
replaced by u' and v'. If the scale factor e - 2 A t ( M +") j n f |D.34j) is expanded as 
e -2fi(tt*+u*) _ 2jj,{u' — u* + v ' — f*)] then the integrals inside the square root in 
(ID. 411) can be broken into two parts: 



(± Lo f e -^*+v*) ( r [ Vl (u>) _ Vm ] du ' + [ U * + ~ AU [v* - v 2 (u')\ du' 

(D.42) 

and 

/ / fiAdv'du'+ / / piAdv'du' 

J u* — 2-k J V* J u* J u'— 2lT— u*+v * 

(D.43) 

where 

A = -2(u' -<)-{F u , -«(«'- |tt)} + (?/ -v*)d u ,Q(u') 

-8{- \n(W + 2u' - 2tt) + ln(2«') + K(u') cos 2u' + £(«') sin 2u'}. 

(D.44) 

Now, it is easily seen that the contribution of the term —2{v! — + v' — i>*) 
in A to the integrals in ( 1D.43I) is zero, and since the remaining terms in A are all 
half-range periodic functions of v! it turns out that 

/•u' + 2TT — U t +Vt, pU t +2-K PVt, 

/ fxAdv' du' + / / fxAdv' du' 

u f ~2irJv* J u t J u 1 — 2-7T— u t +v t 

= 2irfx / (u' — u* — -n)d u iG{u') du — An[i I {F u — k(u — du 

Ju* JO 

-32tt// [ {- \n(W + 2u - 2tt) + ln(2tt') + K(u) cos 2u + £(V) sin 2u'} du. 
Jo 

(D.45) 

Furthermore, the contribution of the second integral here is zero, as -F(vr) = -F(O) 
and u—\t[ integrates to zero. 

Turning to the integrals in (ID. 421) . the terms in V\(u') — v* and t>* — t>2(V) of 
order /i° contribute 4-7T 2 to the integrals as expected, while the terms proportional 
to Au and Av integrate to give zero to order /1. Also, since Q is half-range periodic 
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in u', and n!F(u'] u*) — pJ-{v! — 2n; u*) = —Av from (ID. 33ft and f]D.38[) . it turns 
out that 

u t put+2ir— Ait fUi,+2iT 

[vi(u')—vJdu'+ [u*— v 2 (u')] du = 47r 2 +27r/i / Q(u) du'—2irAv. 



It* — 27T + Atl 



Combining flD.45j) and (ID. 461) gives: 



(D.46) 



•Vl(u') />tt*+27T— Alt pV* 

^ 2 dv'du'+ / / ^dv'du' 

, — 2tt+Aix J J u, J V2(u') 

= 47r 2 + 47r 2 ^K) +27r / u[167r-8Si(27r)] 

-32nfi / {- \n{W + 2u' - 2tt) + ln(2u') + K(u') cos 2u' + £(«') sin 2u'} du' 
Jo 

(D.47) 

where the value of the integral in the final term is 27rln27r + Si(27r) — 3tt. 
So, to first order in /z, the invariant length of the loop is: 

LiKO = 2(^L )e-^ (u * + ^ ) {4vr 2 + 47rV[^ + 4( 7E -Ci^ + lniy,)+Ci] 
-167T^[-87r + 4vr ln2vr + 3 Si(27r)]} 1/2 
= L e -M«.+«.){i + ^[1^ + 2 ( 7e _ Ci W* + In W*) + Cx] 

-2/i[-8 + 41n27r + 37r- 1 Si(27r)]} (D.48) 

for G (0, 7r). For values of outside this range, the term multiplying the 
scale factor e ~ K ^( u *+ v *) should be replaced by its half-range periodic extension. 
Furthermore the w^-dependent part of this term, \ku* + 2(7 E — CiW 7 * + lnW*), 
tends to 2(7 E — Ci27r + ln27r) as — > + , and to as — > tt~. In view of 
the value (16.601) deduced for n in Section 6.3, the two limits are the same, and so 
L\ is continuous for all values of and v*. 
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LIST OF FIGURE CAPTIONS: 

Figure 1: y-z projection of the Allen-Casper-Ottewill loop. 

Figure 2: Schematic representation of the evaporation of the ACO loop. 

Figure 3: Schematic picture of the relationship between the field point [t, x] 
and its image point [i, x] . 

Figure 4a: Variation of /i tt /(4//) with Z along the z-axis. 

Figure 4b: Variation of h tx / (4/i) with Z along the z-axis. 

Figure 4c: Variation of h ty / (4/i) with Z along the z-axis. 

Figure 4d: Variation of h zx / (4/i) with Z along the z-axis. 

Figure 4e: Variation of h zy / (4/i) with Z along the z-axis. 

Figure 5a: The correction G as a function of u if C\ — 9. 514 3. 

Figure 5b: The combination F-t-sgn(ii) A as a function of u if C2 ~ 15. 403. 

Figure 5c: The "natural" correction function A(u) = —^-(\u\ — 

Figure 5d: The correction F as a function of u if A has its natural form. 

Figure 6: The Kibble- Turok sphere of the ACO loop and two possible types 
of perturbation. 
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